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Minimum Fluidizing Velocities 


C. G. SINCLAIR? 


for Coal in Air’ 


This paper reviews the methods which have been proposed for 
evaluating the lower limit of velocity for which a fluidized bed of solids 
will be stable. Experimental tests for determining the minimum fluidizing 
velocity have been obtained for nine narrowly sized coal samples ranging 
from 0.0099 to 0.224 in. in diameter, and for one widely sized sample 
less than a No. 10 Tyler sieve using dry air at atmospheric conditions in 
a 5! in. diameter bed. The results are compared with those predicted 
using the correlations of other workers. 


A a preliminary step in the design 


of a fluidized bed reactor it is 
essential that the limits of the gas 
velocity within which the bed will be 
stable should first be determined. The 
study described in this paper was 
concerned with the lower limit of 
velocity, which is called the critical 
or minimum fluidizing velocity. The 
question arose in connection with a 
project for the carbonization of coal 
in a fluidized bed, when it was 
necessary to determine the minimum 
fluidizing velocity for particles larger 
than those generally studied. 


Definition of the critical fluidizing 
velocity 

Leva et al.(') have defined the 
critical fluidizing velocity as the 
velocity which is obtained when the 
pressure drop is set equal to the 
weight of the bed in the equation 
they propose for corre'ating velocity 
and pressure drop across a packed 
bed. Since their original equation con- 
tains a function of the void fraction 
in the bed this quantity has to be 
determined experimentally at the 
critical point. 
Gilliland and 


Lewis, Bauer (?) 


plotted a curve connecting the system 
properties and bed voidage with par- 
ticle Reynolds number for both the 
fixed and fluidized state. They state 
that the two curves can be solved 
simultaneously for the gas velocity 
and void fraction at the transition 
between the fixed and fluidized state. 
However, this gives imaginary values 
for the void fraction at certain values 
of the Reynolds number, as is shown 
later. 

Baerg, Klassen and Gishler in their 
studies of heat transfer in fluidized 
beds ‘*) noted a sharp discontinuity in 
the plot of heat transfer coefficient to 
the walls of the retaining vessel versus 
superficial gas velocity. They assumed 
that this gas velocity ‘corresponded to 
the transition between the fixed and 
fluidized conditions and were able to 
correlate it with system properties. 

Agarwal and Storrow ‘*) deter- 
mined the critical fluidizing velocity 
by visual observation of the change 
from fixed to fluidized bed condition. 

A graph of pressure drop versus 
velocity forms the basis of two defi- 
nitions of critical fluidizing velocity. 
It is found from this plot that during 
the fixed bed condition the pressure 
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and D. B. ROBINSON: 


drop increases regularly with increase 
in gas velocity. Near the onset of 
fluidization, the curve falls off, and 
during fluidization the pressure is sub- 
stantially constant. If the velocity is 
reduced slowly, the curve will be 
retracted along the same path in the 
fluidized region but will fall off at a 
higher velocity. Thus, in the fixed 
bed region for any given velocity the 
pressure drop will be lower when the 
velocity is being reduced than when 
it is being increased. These features 
are illustrated in Figure 1. The lower 
curve is obtained at maximum voidage 
in the bed and if there is any jarring 
of the equipment or if the velocity is 
reduced suddenly the points will be 
above this line. The slope of the in- 
creasing line depends upon the initial 
packing of the bed and consequently 
if the horizontal and increasing line 
are extended until they intersect, the 
value of V, will depend on the initial 
packing of the bed. Jolley and 
Stanton (>) use this point to define the 
critical fluidizing velocity provided 
the bed is mechanically shaken. How- 
ever, the point Vie of intersection of 
the decreasing and horizontal curves 
is reproducible since the bed is always 
settled in a constant manner and at a 
maximum voidage, and represents a 
velocity above which the bed will 
alwavs be fluidized. This definition, 
independently proposed by Van 
Heerden et al ‘®) and Miller and 
L ogwinuk (7) is superior to the others 
since it is reproducible and since it is 
directly obtained from experimental 
data. It is the one adopted in the 
present w ork. 

Difficulties arise with this definition 
when dealing with wide size range 
mixtures since the pressure drop 
versus velocity curve obtained has the 
shape shown in Figure 2, and there 
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Figure 2— Generalized pressure drop-velocity relationship for widely 


fractions. 


is some doubt as to the point of inter- 
section of the lines in the range V, 
to Vir. It was noticed that while the 
velocity was being increased there 
was still some segregation of larger 
particles at velocities greater than Vn 
though after thorough fluidization 
there was no segregation until the 
velocity dropped below V,,;. Accord- 
ingly Vie was taken as the critical 
fluidizing velocity in this study. 


Review of correlations for the 
critical fluidizing velocity 

The equations connecting pressure 
drop and gas velocity in fixed beds 
can be modified to apply to fluidized 
beds by substituting the identity 


AP 
-- 


By a correlation of the bed voids at 
the critical fluidizing velocity against 
the particle diameter, it is possible to 
solve for the critical fluidizing veloc- 
ity. This was the procedure adopted 
by Leva et al“). It results in the 
following equation: 


€mt) (Ps — pr)..- (1) 


g Dp pi (Ps — pr) €8 
OF As" (1 = €)2 n 


G? = 


Plots of f and n versus Re, are 
given in the paper. By setting « = emg, 
the fraction voids at the critical fluid- 
izing velocity G = Grn = Vingpp May 
be obtained. Plots of émt Versus D,, are 
given by Leva et al. for sand (s)' ‘and 
for anthracite coal (®), 

Van Heerden, Nobel, and Van 
Kreveln ‘®) studied fluidization only 
in the laminar region where Re, < “10. 
For fixed beds the pressure dome can 
be evaluated by an equation of the 
type 


AP fg (€) B.2p; V? 
L D, 


(3) 


Gilliland and 


¢@ 


For Re,<10, Lewis, 
Bauer (2) have found that f Re. 
and for the fluidized bed one can set 
AP 

L 

From this it follows that 


= £ Pom 


Pt Pom & D,? = Vt Pr Dp 


Rems es oe 
154 Me g (€me) p KM 


This equation relates the critical 
fluidizing velocity V,,~¢ to the fraction 


voids at the critical,e,,,. Van Heerden 
found that 6 as determined using this 
equation varied with D, for sub- 
stances with micro- porosity such as 
coke. However, both his data and that 
of Lewis, Gilliland and Bauer ‘) indi- 
cate that spheres have a definite 
porosity of 0.406 at the critical point 
and this porosity is independent of 
D,. It seemed logical to postulate, 
therefore, that with substances havi ing 
micro-porosity such as coke and for 
irregular particles, the effective void- 
age for gas flow would also be 
independent of D,. As this effective 
voidage cannot be measured, Van 
Heerden suggests that g(en¢-) and 6 
be combined to give a new shape 
factor B. This can now be determined 
experimentally by comparing the crit- 
ical fluidizing velocity of the irregular 
particles w ith that for spheres having 
the same nominal diameter. By setting 
B equal to 1 for spheres, the followi ing 
equation was obtained from the ex- 

perimental results: 


0.00123 ¢ 


———— Pom Pi Dy?... . .(5) 
B gy? 


Remi = 


Values of B evaluated using this equa- 
tion show much less spread than the 
corresponding § values. Though this 
equation holds for particles of uni- 
form sizes they found that for mix- 
tures of wide range the shape factor 
B is sti'l not constant with the arith- 
metic mean diameter of the particles. 
If the equivalent diameter of the par- 
ticles is redefined as the diameter of 
those spheres with the same number 
per unit volume of the bed as the 
particles, both measured at maximum 
porosity, then the discrepancy is 
eliminated. 


Agarwal and Storrow ‘*) found that 
the fluidizing gas velocity is propor- 
tional to the terminal velocity of the 
particles calculated from a mean 
diameter (equal to six over the par- 
ticle surface area) and proportional to 
the bed voidage raised to the sixth 
power. They show plots of the 
critical bed voidage against particle 
diameter for the materials they 
worked with, thus the critical fluidiz- 
ing velocity can be determined. They 
found that their experimental results 
did not check with values calculated 
from Leva’s correlation. 

Miller and Logwinuk 7)  ap- 
proached _ the problem from the 
dimensional analysis point of view, 
which resulted in the following equa- 
tion involving particle and gas pro- 
perties raised to suitable powers: 


0.00125 
Mu 


D,? (ps — pi )0-9 pelt g. (6) 


The solids used in their investigations 
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were silicon carbide, aluminum oxide, 
silica gel and sand. The equation is 
valid in the laminar flow region. 
Lewis, Gilliland and Bauer (2) in 
the course of their study of fluidiza- 
tion in the laminar region obtained 
two plots of friction factor (which 
includes system properties and gas 
velocity) for the fluidized bed and 
the fixed bed, against particle Rey- 
no'ds number. They state that these 
two curves can be solved simultan- 
eously for the critical fluidizing 
velocity and fraction voids. These 
friction factors are respectively, for 


the fixed bed 
g Dp (AP/L) 
2 V2 pt 


3 é& 


"(1 — €) 


B (i — 
and for the fluidized bed 


4g Dy (me — m)ee™ 5 ues 


3 V? p: 


For a given system of particles and 
gas, the ratio of the fluidized to fixed 
bed friction factors becomes 


fr 8 (ps — pi) 
fg 3 (AP) 
L 
and assuming that the relation 


AP 
* (Ps — pr) (1 — €) 
is valid at the critical point the ratio 


further reduces to: 


fr 

fs 3 
Hence by determining from the ex- 
perimental results the va-ue of this 
ratio at different Reynold’s numbers, 
it should be possible’ to obtain a plot 
of the fraction voids at the critical 
point versus Reynold’s number for 
any system of ‘spheres and a gas. 
These’ computations have been per- 
formed and are tabulated in Table 1. 
It appears that this method cannot be 
used for determining the fraction 
voids or the gas velocity at the critical 
point. 

Wilhe'm and Kwauk (!°) have pub- 
lished their work on a wide range of 
fluidized systems in a most useful 
form. Their general correlation is a 
plot of the two dimensionless groups. 


it? 
g ce 


Kap = op — 


€!.65 (1 a €)? 


eh = -@) 


Dp? pr & (ps = pr) 


K => 
Ap 2 a? 


against the partic'e Reynold’s number 
Re,. As a fluid stream is first passed 
up through a bed of particles, condi- 


TABLE 1 


EVALUATION OF € AT THE CRITICAL POINT OF A FLUIDIZED BED 
(after Lewis, Gilliland and Bauer) 


tions are equivalent to a fixed bed and 
one curve represents the general cor- 
relation between KAp and Re,. At 
the critical point according to equa- 
= (ps — pr) (1 — e) and 

Rea AP 


a ae : so that | = 


(ps — pr) and KAp = Kap. They 
found that this is in fact the case and 
the individual curves break off hori- 


tion (1) 


—4 + 8 Tyler 


TOT 


+ 14 Tyler 


—50 + 65 Tyler 


Figure 3—Microphotographs of coal 
samples compared to scale in milli- 
metres. 
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zontally from the first curve and are 
terminated on or before a second 
curve, which represents the relation- 
ship between KA p and Re, for free 
falling particles. 

Pinchbeck and Popper ‘!!) studied 
the fluidization of spherical and non- 
spherical particles using nitrogen. By 
considering the total force holding 
a particle in suspension to be the 
sum of viscous resistance and fluid 
impact they show that the ratio of 
the value of the particle Reynold’s 
number at the terminal velocity to the 
value at the critical velocity is a func- 
tion of the group 

py? 
g (ps — pr) px Dp? 
It is not stated how the critical 
ve ocity was determined. 

Zenz (12) presents curves comparing 
the work of previous investigators on 
the basis of a drag coefficient using 
porosity as a parameter which has a 
range from about 0.4 at minimum 
fluidization to 1.0 at the Stokes 
velocity for the particles in the bed. 
His calculations show that the distri- 
bution of sizes in a bed of mixed 
sizes is an important consideration in 
choosing the effective particle diam- 
eter for calculating bed characteristics 
for fluidization. 

Zaki \'3) reports the results of ex- 
tensive tests on the fluidization of a 
variety of particles using water, and 
he has evaluated the variation of the 
porosity function with Reynold’s 
number and the ratio of particle to 
tube diameter. 

Fina ly, Baerg, Klasen and Gishler 
(2) propose the empirical equation. 


Gm = 1.3 x 108 Dy pp: 


determined from their heat transfer 
data. 
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Figure 4—Pressure drop in a packed and fluidized bed of 
—4 + 8 crushed coal. 
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Figure 6—Pressure drop in a packed and fluidized bed of 
—48 + 65 crushed coal. 
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Figure 5—Pressure drop in a packed and fluidized bed of 


—10 + 14 crushed coal. 
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Air Flow Rate -Standard Cubic Feet Per Second 
Figure 7—Pressure drop in a packed and fluidized bed of 
—10 crushed coal. 
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TABLE 2 
PROPERTIES OF COAL SAMPLES TESTED 




















Critical Fluidizing Critical Reynolds No. Density 
Tyler Arithmetic | Geometric Velocity Equivalent based on lb/ft. 
Sample sieve Mean Mean = diam. 
limits diam. | diam, inches 
in. in. Ib/hr/ft.? ft/sec. A.M. Equiv. bulk | true 
| | diam. die im. 
 * ~3 +4 | 0.22 | 0.220 | 1,437.0 | 5.75 615.0 44.1 | 87.6 
2 -4 +8 | 0.139 | 0.131 | 1,005.0 | 4.28 0.129 267.0 | 249.0 | 43.5 | 87.6 
3 —8 +410 | 0.079 .0776 718.0 2.92 0.0858 108.3 117.7 42.0 87.6 
| 
+ -10 +14 | 0.056 0546 | 526.0 2.10 0.0628 56.3 63.1 41.1 87.6 
5 —14 +20 | 0.0394 | .0388 | 348.5 1.425 0.0428 26.3 28.5 | 40.2 87.6 
6 —20 +28 | 0.028 | 0276 | 210.2 0.841 0.0301 11.26 12.1 | 39.5 87.6 
| 
7 —28 +35 | 0.0198 | 0195 | 103.7 0.417 | 0.0211 3.93 4.19 | 38.8 87.6 
8 —35 +48 0.014 £0135 50.1 0.203 0.0136 1.34 4:31 | 38.0 87.6 
9 —48 +65 0.0099 .00975 | 26.7 0.108 0.0100 0.505 32-2 87.6 
- i = sci 
Experimental and an average reading was deter- TABLE 3 


The equipment used for determin- 
ing the critical fluidizing velocities of 
the coal was of standard design. Dried 
air, which was metered by an orifice 
plate, was passed through the fluidiza- 
tion equipment. This consisted of an 
18° metal cone, a porous bronze sup- 
porting screen for the bed of particles, 
and a 54 inch I.D. Lucite tube. The 
pressure drop across the bed was 
measured from two pressure taps, one 
immediately upstream of the porous 
bronze screen and one above the 
fluidized bed. The pressure drop 
versus air flow rate curve for the 
empty tube was determined prior to 
making the fluidization tests. 

A weighed amount of the coal 
sample was charged to the fluidization 
tube, and the correct orifice plate to 
cover the range of air rates antici- 
pated was put in position. The equip- 
ment was closed up and the coal 
sample was gently fluidized to even 
out the surface. By manipulation of 
the control valves the air rate was in- 
creased in steps, allowing sufficient 
time for equilibrium at each step. 
After each stepwise change in air 
rates, the temperature, static pres- 
sures, differential pressures, bed 
height and bed appearance were 
noted. For each run the air rate was 
increased in steps up to and beyond 
fluidization, then reduced in steps to 
zero. 

After fluidization the bed height 
readings were less accurate and when 
the bed began to slug an approximate 
average was all that could be read. 
Similarly after fluidization the mano- 
meter levels would oscil‘ate slightly 


mined by eye from the maximum and 
minimum values. 


Experimental results 


The fluidization characteristics and 
critical fluidizing velocity for beds 
of coal im air were obtained experi- 
mentally for nine different coal sizes 
between No. 3 and No. 65 Tyler 
mesh. The samples were prepared by 
crushing a sub-bituminous B_ coal 
from Edmonton, Alberta, in a ham- 
mer mill and sizing between Tyler 
sieves. 


Samples of the coal particles used 
for the fluidization tests are shown in 
Figure 3. The comparative scale is 
divided in millimeters. The physical 
properties of these particles are 
reported in Table 2. The equiva’ent 
diameter was calculated by the 
method proposed by Van Heerden (°). 


Typical examples of the experi- 
mentally determined curves are shown 
in Figures 4, 5, 6 and 7. Figure 4 
indicates the data obtained for the 
—4 + 8 fraction, and Figure 5 the data 
for the —10 + 14 fraction in beds 
containing 3 lb. of coal. The data for 
the increasing flow rates in Figure 5 
were obtained for the coal as it was 
poured into the container, in order to 
indicate the pronounced effect of the 
initial state of the bed on pressure 
drop. Figure 6 shows the data for the 
—48 + 65 fraction in a bed contain- 
ing 5 lb. of coal. Figure 7 represents 
the data for a sample of mixed sizes 
less than Tyler mesh 10. The screen 
analysis of this sample is given in 
Table 3. It will be noted that the 
characteristic hump _ reported by 
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SCREEN ANALYSIS OF WIDELY SIZED 
CoaAL FRACTION 


Tyler Screen Percentage 
0 + 0.4 
—14 +20 11.1 
—20 +28 18.8 
—28 +35 17.4 
—35 +48 14.4 
—48 +65 $.3 
—65 +4100 ms 
—100 +150 6.4 
—150 +200 3.8 
—200 10.7 


many workers is absent in the larger 
sizes and in the mixed sample. The 
results of all tests, including the ex- 
perimental data on which the curves 
in these figures are based, are shown 
in Table 4. 


Calculated results 


The experimentally determined 
critical fluidizing velocities were 
compared to the calculated velocities 
using the correlations of Van 
Heerden (6), Leva et al “), Miller and 
Logwinuk (7), and Baerg, Klassen, and 
Gishler (3), Of these four correlations, 
those of Miler and Logwinuk, Equa- 
tion 6, and Baerg, Klassen, and 
Gish!er, Equation 7, present no diffi- 
culties. 


55 





Run 1 Sample 4 | Run 2 Sample 4 


Wt. Charged Wt. Charged Wt. 
5 Ibs. 3 Ibs. f: 
T = 530°R, T = §29°R, T = 
P = 703 mns. P = 706 mns. P = 
Air Pres- Air Pres- Air 
Flow sure Flow sure Flow 
Rate Drop Rate Drop Rate 
0.068 2.25 0.088 2.60 0.107 
0.128 2.95 0.132 5.00 0.159 
0.139 5.30 0.178 7.05 0.198 
0.206 9.50 0.215 8.90 0.224 
0.232 | 11.20 0.225 8.45 0.250 
0.255 2.75 0.238 8.75 0.279 
0.266 13.40 0.248 8.70 0.292 
0.274 | 13.70 0.258 8.85 0.304 
0.282 | 13.95 0.269 8.65 0.325 
0.289 | 14.10 0.279 8.40 0.350 
0.296 | 14.05 0.289 8.30 0.380 
0.307 |; 14.15 0.297 8.35 0.376 
0.329 14.5 0.314 8.45 0.309 
0.359 14.6 0.332 8.45 0.277 
0.319 14.25 0.356 65 0.229 
0.292 | 13.65 0.374 8.70 0.187 
0.261 11.50 0.325 8.50 0.139 
0.208 8.75 0.307 8.25 0.125 
0.144 5.50 0.300 8.10 0.188 
0.160 6.15 0.288 7.70 0.224 
0.225 9.85 0.276 | 7.30 0.259 
0.264 12.45 0.255 6.55 0.282 
0.278 | 13.40 0.228 3.33 0.302 
0.286 | 14.00 0.189 4.50 0.328 
0.290 | 0.134 2.75 0.363 
0.294 14.55 0.092 | 1.70 0.397 
0.301 | 14.50 0.146 | 3.10 0.112 
0.306 | 14.45 0.188 4.60 0.171 
0.312 | 14.25 0.236 6.20 0.216 
0.321 | 14.35 0.265 7.40 0.250 
0.346 | 14.60 | 0.277 | 7.90 | 0.270 
| 0.296 8.35 0.283 

| 0.298 | 8.55 | 0.297 





The bulk density of the solid p, 
was obtained by pouring the solid 
under control'ed conditions into a 
measuring cylinder and finding the 
weight of a known volume. The 
density of the solid ps Was determined 
pyknometrically using pure toluene 
and was found to be 87.6 lb/ft® for 
all particle sizes. The density of the 
fluidizing air was calculated using the 
ideal gas law knowing the temper- 
ature and pressure. The particle diam- 
eter was taken as the arithmetic mean 
of the sieve openings bracketing the 
sample for Baerg, Klassen and 
Gishler’s equation, and the geometric 
mean for Miller and Logwinuk’s 
equation. The difference between the 
arithmetic mean and the geometric 
mean is less than 1.5% in all cases. 


In order to determine the shape 
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Run 3 Sample 4 


TABLE 4 
EXPERIMENTAL FLOW RATES AND PRESSURE DROPS 


Run 4 Sample 5 


Run 5 Sample 5 


Run 6 Sample 6 


Run 7 Sample 7 


Charged Wt. Charged Wt. Charged Wt. Charged Wt. Charged 
5 Ibs. 5 lbs. 3 Ibs. 5 lbs. 5 lbs. 
540°R, T = 541°R, T = 534R, T = 537°R, T = 539°R, 
704 mns. P = 703 mns. > = 711 mns. P = 713 mns. P = 711 mns. 
Pres- Air Pres- Air Pres- Air Pres- Air Pres- 
sure Flow sure Flow sure Flow sure Flow sure 
Drop Rate Drop Rate Drop Rate Dr op Rate Drop 
2.00 0.037 2.40 0.099 4.05 0.0308 3.50 0.0169 | 2.85 
2.25 0.074 4.80 0.139 5.85 0.0527 6.28 0.0239 | 5.05 
2.30 0.099 6.95 0.166 7.45 0.0700 8.65 0.0338 | 8.09 
2.65 0.127 9.30 0.194 8.45 0.0836 | 10.55 0.0441 | 11.06 
2.95 0.151 11.70 0.223 8.65 0.0992 | 12.80 0.0541 | 13.42 
3.50 0.171 13.65 0.246 8.85 0.115 14.07 0.0671 | 14.67 
3.60 0.187 | 14.50 0.267 8.95 0.132 14.48 0.0831 | 15.00 
3.65 0.201 14.50 0.288 8.95 0.148 14.60 0.111 | 15.49 
3.70 0.212 14.60 0.307 8.90 0.172 14.58 0.216: | 15,81 
eb 0.246 14.70 0.291 9.00 0.202 14.60 0.144 16.12 
3.95 O,272 | 14.75 0.272 9.00 0.188 14.64 0.125 15.72 
3.65 0.314 | 14.80 0.247 | 8.95 0.164 | 14.58 0.111 | 15.44 
3.90 0.291 14.80 0.224 8.85 0.142 14.54 0.0950 | $5.55 
3.45 0.253 14.75 0.204 8.60 0.125 14.32 0.0739 | 14.65 
2.65 0.219 14.60 0.178 7.50 0.113 13.07 0.0638 | 13.88 
2.00 | 0.203 | 14.40 | 0.145 | 5.85 | 0.0875 | 10.09 | 0.0521 | 11.33 
1.40 | 0.162 | 11.70 | 0.119 | 4.75 | 0.0706 | 7.80 | 0.0421} 8.52 
1.20 0.107 7.30 0.088 3.30 0.0519 | 5.53 0.0329 6.04 
2.10 0.081 5.30 0.128 §.25 0.0351 | 3.29 0.0262 aaa 
2.80 0.108 7.38 0.173 7.45 | 0.0339 | 3.85 0.0267 | 6.17 
5.25 0.134 9.60 0.199 8.40 0.0579 6.80 0.0439 11.80 
| 
3.50 0.153 | 11.50 0.211 8.40 0.0823 | 10.05 0.0488 | 14.44 
3.65 0.174 | 13.40 0.248 8.90 0.110 | 13.85 0.0718 | 14.70 
75 0.187 | 14.35 0.274 | 8.90 0.143 | 14.58 0.0971 | 15.30 
3.80 | 0.190 | 14.60 | 0.179 | 14.65 | 0.128 | 15.85 
3.95 | 0.201 | 14.30 
1.85 0.210 14.55 | 
3.20 | 0.257 | 14.75 | 
3.90 | 0.295 | 14.85 | 
3.80 | 
3.75 | | 
3.80 





factor ) for Equation 2 of Leva et al, 
the coal was compared with photos 
of an antricite coal which Leva (9) 
quotes as having a shape factor of 1.6. 
The two materials appeared to be of 
simi'ar shape and the value 1.6 was 
therefore used. On this basis it was 
assumed that at the flow rate G,,, the 
minimum fluid voidage for a bed 
would be the same as the experi- 
mentally determined values quoted in 
this paper. Values of f and n for 
various ies of Re, were obtained 
from the plots given. 

The difficulty in using the equation 
proposed by Van Heerden et al ‘®) 
lies in estimating py», the bed density 
at the critical point. However, it is 
stated in their paper that the differ- 
ence between the bed density at the 
critical point and that determined by 


filling a graduated cylinder with the 
loosely packed particles is relatively 
small, and recommend the use of a 
value obtained this way. The value 
of the equivalent diameter De may 
be determined from the bed density 
and the number of particles in a given 
weight of material. 


The values of the critical gas 
velocity calculated from these various 
correlations are compared graphically 
to the experimental values in Figure 
8. The correlation of Baerg, Klassen 
and Gishler obviously does not apply 
to this system since the calculated 
velocities range from 30-70% below 
the experimental values. The veloci- 
ties calculated from the correlation of 
Leva et al range from 7 to 30% below 
the experimentally determined values. 
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Run 9 Sample 2 


Rin 8 Sample 3 
Wt. Charged 


Vt. Charged 


EXPERIMENTAL FLOW RATES AND PRESSURE DROPS 


TABLE 4—Continued 


Run 10 Sample 8 


Run 11 Sample 9 


Wt. Charged Wt. Charged 








3 Ibs. 3 lbs. 5 lbs. 5 Ibs. 

lr = 543°R, T = 538°R, T = 536°R, T = 534°R, 

= 709 mns. P = 702 mns. P = 702 mns. P = 702 mns. 
si eae areas he enc Eee coa " 

\ir Pres- Air Pres- Air Pres- Air Pres- 
Fiow sure Flow sure Flow sure Flow sure 
Rate Drop Rate Drop Rate Drop Rate Drop 
0.127 | 1.50 0.094 0.40 0. 00517 2.49 0.00331 2.90 
0.226 ato 0.157 1.30 0.00852 4.98 0.00496 5.39 
0.267 | 5.20 0.207 2.00 0.0134 8.06 0.00670 | 8.08 
0.314 6.45 0. 263 2.80 0.0185 10.85 0.00878 | 11.03 
0.361 7.85 0.329 4.05 0.0239 13.72 0.0112 13.82 
0.422 8.65 0.404 5.30 0.0281 14.56 0.0131 14.46 
0.476 8.90 0.453 6.40 0.0326 14.30 0.0153 13.76 
0.535 | 8.75 0.521 7.75 0.0386 14.42 0.0185 13.85 
0.600 8.85 0.561 8.50 0.0465 | 14.60 0.0215 13.99 

| | 
0.482 8.95 0.592 | 8.70 0.0549 14.62 0.0255 |} 14.17 
0.455 9.00 0.649 | 8.75 0.0696 14.71 0.0294 14.30 
0). 387 7.95 0.724 8.60 0.0800 14.72 0.0348 | 14.45 
0.360 7.05 0.758 8.80 0.0654 } 14.73 0.0310 14.45 
0.330 6.15 0.686 8.75 0.0574 |} 14.70 0.0286 14.36 
0.295 $25 0.627 8.60 0.0513 | 14.68 0.0256 14.32 
0.252 4.10 0.551 | 7.80 0.0434 16.61 0.0236 14.33 
0.189 2.50 0.522 7.10 0.0351 14.44 0.0217 14.24 
0.624 4.95 0.469 6.20 0.0294 13.55 0.0197 | 14.04 
0.345 y ee @.431 | 5.45 0.0260 | 22i@e 0.0179 13.85 
0.392 8.55 0.393 | 4.70 0.0210 | 10.39 0.0149 13.26 
0.414 8.70 0.339 | o.45 0.0132 | 6.61 0.0131 12.36 
0.494 8.95 0.256. | 2.35 0.00875 | 4.58 0.0109 | 10.92 
0.335 | 3.95 0.0152 | 8.15 0.00855 | 8.88 
0.480 | 6.80 0.0221 11.83 0.00539 5.19 
0.676 8.75 0.0275 } 14.27 0.00355 | 2.00 

0.738 8.60 0.0323 } 14.15 | 

0.0438 14.41 | 

| 

| 0.0571 | 14.66 


The two quantities used in the cor- 
re'ating equation which are in doubt 
are the shape factor ) and the porosity 
at the critical point ¢,,;- In order to 
evaluate the influence these quantities 
have on the discrepancies, the values 
of 4 and ene necessary to equalize the 
experimental and calculated values 
were separately determined. This was 
done by calculating the value of A 
necessary to make the predicted and 
experimental curves coincide when 
using Leva’s values for ¢~. The aver- 
age value of the shape factor so deter- 
mined was 1.20 and it showed a mini- 
mum of 1.01 where the diameter was 
0.056 in. Similarly, values of ¢.,, were 
obtained for a fixed value of A of 1.60 
and were found to range from 0.51 to 
0.60. It is not thought that the original 
estimates were this much in error, 


rather it is thought that these dis- 
crepancies rising to a maximum in 
the intermediate flow range are due 
to the inaccuracies inherent in apply- 
ing a friction factor determined from 
fixed bed measurements to fluidized 
bed conditions. 


The correlation of Van Heerden et 
al (6) gives critical fluidizing velocities 
for the two sma'lest size of particles 
20% too low and for the two larger 
sizes 7% too low. This may be due 
to inaccuracies in determining the 
equivalent diameter for the smaller 
particle sizes. It will be noted that for 


the two smallest sizes there is no 
significant difference between the 
arithmetic mean diameter and the 


equiva'ent diameter of Van Heerden. 
The average value of the shape factor 
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Run 12 Mixed Sample 


Run 13 Sample 1 


Wt. Charged Wt. Charged 





3.8lbs. 3 Ibs. 
T = 542°R, T = 532°R, 

P = 707 mns. P = 704 mns. 

Air Pres- Air Pres- 

Flow sure Flow sure 

Rate _Drop Rate Drop 
0.00583 1.20 0 253 0.77 
0.00729 3.58 0.440 2.10 
0.0115 6.06 0.573 3.12 
0.0163 7.95 0.684 4.05 
0.0230 10.28 0.860 5.85 
0.0358 10.89 0.997 5.70 
0.0424 11.26 1.084 5.50 
0.0514 | 11.23 1.153 4.75 
0.0582 11.40 1.250 4.95 
0.0645 | 11.38 1.400 | 4.60 
0.0595 | 11.40 1.533 | 4.30 
0.0533 11.32 1.402 4.50 
0.0467 | 11.25 1.235 5.55 
0.0386 | 10.87 1.139 5.65 
0.0338 | 10.19 1.012 6.20 
0.0272 | 9.12 0.869 6.60 
0.0144 7.14 0.669 4.50 
0.0144 5.66 0.431 2.20 
0.00974 3.67 





B which would cause the experi- 
mental and calculated values to coin- 
cide was determined and found to be 
0.8. 

The correlation of Miller ‘and 
Logwinuk ‘7) gives the most satisfac- 
tory results but its application is 
restricted to the laminar region. The 
calcu'ated values average 9.5% low 
but agree with the experimental if the 
constant 0.00125 in Equation 6 is 
changed to 0.00138. 

The influence of the choice of the 
critical fluidizing velocity is demon- 
strated in Figure 9 by comparing the 
experimental data to the curves of 
Wilhelm and Kwauk “), They ob- 
tained their data as the fluid velocity 
was being increased and hence the 
values do not correspond to the criti- 
cal velocity as defined here. 
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Figure 9—Comparison of results with data of Wilhelm and Kwauk (!°), 
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Conclusions 


The results of this work have 
shown that differences exist betwcen 
experimental and ca ‘culated critical 
fluidizing velocities when app: ying 
the correlations proposed by several 
workers to the case of fluidizing ir- 
regular shaped particles having micro- 
porosity such as coal. They have also 
pointed to the lack of av ailable cor- 
relations for the turbu'ent region 
which is encountered when fluidizing 
particles somewhat larger than usual. 
It is most important to choose a mini- 
mum fluidizing velocity that is repro- 
ducible when correlating the experi- 
mental results. 


Nomenclature 

B = Factor allowing for shape and bed 
por« sit y 

D, = Particle diameter, ft. 

De = Equivalent particle diameter, ft. 

f = Friction factor 

fs = Friction factor for fixed bed 

fr = Friction factor for fluidized bed 

g = Acceleration of gravity, ft/sec.? 

g (€) = Porosity function 

G = Mass velocity of fluidizing gas, 
lb/ft? sec. 

L = Height of bed, ft. 

Lo = Height of completely compacted 
bed, ft 

n = Power used in Leva (1) correlation 

AP = Pressure drop, lb/ft.? 

Re, = Reynold’s number based on par- 
ticle diameter 

V = Superficial gas velocity, ft/sec. 

B = Factor allowing for particle shape 

€ = Bed porosity 

ps Density of the solid, lb/ft. 

pi = Density of the fluidizing gas, 
lb/ft. 

Pov = Bulk density of the bed, lb/ft.’ 

Ppm = Bulk density of the bed at point 
of minimum fluidization, Ib/ft.3 

= Particle shape factor 

be = Fluid viscosity, lb/ft. sec. 

mf = Subscript denoting condition at 


point of minimum fluidization 
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Case Study of Drives for 


Petrochemical Compressors’ 


BRUCE BIRRELL? 


The characteristics of drivers for large petroleum refinery com- 
pressors are discussed. Gas turbines, steam turbines and electric motors 
are compared for six different applications. If cheap refinery gas is 
available the gas turbine has an advantage but otherwise electric drivers 


are favored. 


<? petrochemical project at the 
Sarnia refinery of Imperial Oil 
Limited was undertaken to supply 
the Canadian market with basic raw 
materials for the manufacture of the 
numerous petrochemical products on 
sale today. The two most important 
materials which will be produced are 
ethylene and butadiene. Ethylene is 
used in the manufacture of nylon 
intermediates, synthetic rubber, latex 
base paints, etc. 

Imperial’s petrochemical plant con- 
sists of a steam cracker, a fraction- 
ating unit, an ethylene recovery 
plant, a butadiene extraction unit, 
and a clay treating plant. The plant 
is designed to crack and process a 
variety of hydrocarbon feed stocks 
ranging from a heavy gas oil to a 
light naphtha. 

One of the important steps in the 
process is gas compression. This is 
important because the essential pro- 
ducts are produced in the gaseous 
phase. In order to recover and purify 
these hydrocarbons by fractionation 
they have to be compressed from the 
relatively low operating pressure of 
the cracking process to the high 
pressure necessary to liquefy them 
for their fractional separation. 

Sections 1-4 below are an appraisal 
of gas turbines, steam turbines and 
electric motors individually, and as 
they apply to driving large com- 
pressors. 

Section 5 (tables 3 to 7 inclusive) 
is an economic study of a particular 


industrial application to ascertain 
which type or types of driver to 
employ. 


Section 6 gives the conclusions 
reached from the study. 


1. Drivers for large compressors 


The main drivers required are, for 
process gas compression and for re- 
frigeration. The propylene compres- 
sion required 3,200 horsepower and 
the remainder 11,600 horsepower, 
totalling 14,800 horsepower. 


Before considering the main fac- 
tors governing the selection of 
drivers, let us consider some of the 
characteristics of the compressors 
themselves that will affect the choice 
of drivers. 

For the large process gas com- 
pressor either reciprocating or cen- 
trifugal machines can be considered. 
Generally, the larger the machine 
the more suited is the centrifugal 
compressor. These machines, how- 
ever, have their limitations. They are 
not suitable for inlet and exhaust 
volumes less than 1,500 CFM and 
500 CFM, respectively, and a gas 
horsepower of less than 500. Also in 
their simplified form they are not 
normally stable in operation below 


50-60% of full load. (This latter 
limitation can, however, be reduced 
in various ways such as by fitting 


guide vanes). “Also, in the case of 
refinery gas, complicated seal sys- 
tems are required to prevent the gas 
under compression from contaminat- 


1Manuscript received November 20, 1957. 


2Manufacturing Department, Engineering Division, 


Imperial Oil Limited, Sarnia, Ont. 


Contribution from the Engineering Division, Impe rial Oil Limite d, Sarnia, Ont. 
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ing the lubricating oil and the air in 
the compressor room. 

On the other hand the large cen- 
trifugal compressor offers low first 
cost and good mechanical balance, 
although the overall efficiency of a 
centrifugal machine is not as great as 
that of a reciprocating machine. 


A reciprocating compressor is in- 


herently unbalanced and may re- 
quire a foundation five or more 
times its own weight. Also the re- 


ciprocating compressor has extensive 


wearing surfaces, cylinder liners, 
links and cross heads. Maintenance 
can, therefore, be a major factor 


particularly when handling corrosive 
gases. Apart from maintenance con- 
siderations the valves associated with 
the reciprocating compressor may 
cause operational difficulties if the 
gas under compression is corrosive 
and causes them to corrode or stick. 
A large reciprocating compressor is 
more suited to a reciprocating type 
driver. A steam turbine driver would 
require a reduction gear which could 
be subject to excessive wear on ac- 
count of the cyclic torque variations 
of the reciprocating compressor. 


Centrifugal compressors are similar 
to centrifugal pumps and are gen- 
erally designed to develop a head of 
10,000 feet per stage compared with 
400 feet per stage in the case of a 
fluid pump. A centrifugal compres- 
sor operates at speeds of the order 
of 3,000 to 9,000 RPM in order to 
develop the required pressure and 
will generally require several stages. 
This results in a high inertia com- 
pared w ith a centrifugal pump, and 
a higher percentage of the total 
inertia of the compressor and driver 
than in the case of pump and driver. 
This introduces problems in starting, 
in the case of electric 
This study is based 


particularly 
motor drivers. 
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on the use of centrifugal com- 


pressors. 

Since the capacities and pressures 
of most of the services in the case 
under review fall within the lim’ts 
quoted above, centrifugal compres- 
sors were selected in preference to 
reciprocating machines. 

The main factors to be considered 
when selecting types of drivers for 
these compressors are: 

(a) The technical requirements of 
the driven units such as horse- 
power, speed, speed range. In 
the case of large units, further 
technical requirements are iner- 
tia, the magnitude of the start- 
ing current in the case of elec- 
tric drivers, the time taken to 
run up to full speed, and, of 
course, such considerations as 
space and weight. In the case of 
back pressure turbines, the 
ability of the system to absorb 
the exhaust steam must be con- 
sidered. 

(b) Public Relations — Noise, the 
purchase of Canadian and locally 
made equipment, and the em- 
ployment of local labor when- 
ever possible. 

(c) Safety as affecting personnel. 

(d) The availability, and cost, of 
fuel or power. 

(e) The overall economics of op- 
eration — depreciation, fuel op- 
erating and maintenance costs. 

(f) The degree of reliability of the 
driver demanded by the process 
particularly in the case of large 
drivers. This is a function both 
of the drivers and the sources 
of fuel on which they depend, 
and is also closely related to 
maintenance requirements. 

Other factors which may have 
a bearing on the selection of the 
driver are: 

(g) Delivery time. 

(h) Availability of replacement 
spares and factory trained ser- 
vice, in the case of specialized 
drivers such as gas turbines. 

(i) Experience with similar drivers 
at existing refineries in the past. 

Although an attempt has _ been 
made to arrange the above factors in 
priority order, it is obvious that 
priorities will vary with each unit 
considered. For example, although 
overall economy of operation is a 
prime factor, it would not take pre- 
cedence over reliability in the case 
of a driver whose failure would 
jeopardize the operation of a com- 
plete process unit. 


With regard to factors (b) and 
(c), it is only by working with a 
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particular organization that one can 
fully appreciate the time and money 
spent in the interests of safety and 
of public relations. 


A particular project is generally 
reviewed on two broad bases: One, 
assuming mostly electric drives, the 
other, assuming mostly steam drives. 
Steam is usually considered to be the 
more reliable of the two, and is 
therefore preferred for important 
drivers such as key process pumps 
(hot and cold charge pumps), fluid 
cat. cracker blowers and compres- 
sors, boiler feed pumps and fan tur- 
bines, cooling water and fire pump 
turbines. 

Electricity, however, offers a low 
capital cost with easier operation and 
maintenance and has the added ad- 
vantage of not requiring a govern- 
ment certified operator. About the 
end of World War II there was a 
distinct trend toward electric drives, 
mainly because of low capital cost 
and secondly due to improved reli- 
ability of electric power supply. In 
the case of steam generated within 
a refinery the employment of high 
operating temperatures and _pres- 
sures, complicated combustion con- 
trol equipment, and only one or two 
steam generators can lead to more 
steam failures than would at first be 
expected. 

Such process units as fluid cat. 
crackers and powerformers have one 
or two comparatively large drivers 
of from 2,000 to 8,000 horsepower 
each, driving compressors on which 
the operation of the whole process 
depends. Accordingly, reliability be- 
comes a prime factor and for that 
particular driver can outweigh in- 
vestment costs. The requirement of 
such process units is for 11 to 15 
months non-stop operation. 

An important point in obtaining 
long runs without a shut-down in 
the case of steam turbine drivers is 
the purity of the steam, and the 
absence of silica with its turbine 
blade fouling tendencies. The em- 
ployment of steam pressures and 
temperatures of 600 p.s.i.g. and 
750°F., and the absence of conden- 
sate recovery calls for a high degree 
of chemical feed treatment of boiler 
water, or alternatively the installa- 
tion of turbine on-load or off-load 
blade washing facilities. This may re- 
quire the sparing of large turbine 
rotors resulting in a _ considerable 
capital outlay. Whenever time per- 
mits, it is well worthwhile to develop 
simultaneously an overall economic 
study that takes into account con- 
densate recovery, chemical feed 
treatment, operating temperatures 


and pressures, turbine washing facili- 
ties and the sparing of large turbine 
rotors. 


The cost of an unintentional shut- 
down of a large unit is difficult to 
assess. It can depend on the percent- 
age of process unit throughput at the 
time, the amount of product in stor- 
age and the ability of the process 
unit to make good the lost produc- 
tion incurred. A large process unit 
operating at a fraction of its design 
rate may make good a two day shut- 
down at nominal cost; on the other 
hand the lost production may 
amount to a thousand dollars per 
hour, or more. 


Development of the industrial gas 
turbine has resulted in its being in- 
cluded in the preliminary study for 
the petrochemical project. No doubt 
in the not too distant future we shall 
automatically include atomic energy 
in our studies of drivers. One 
thought to be borne in mind, espe- 
cially by the utility engineer with 
an oil company is ‘that where ever 
possible our capital should be in- 
vested in process equipment and the 
selected drivers should use a fuel 
obtainable within our own plant. 
Thus, whenever practicable and 
economical, we should think in terms 
of diesel, gasoline and gas engines, 
the gas turbine and steam generated 
from an oil product, rather than in 
terms of electricity or atomic power. 


2. The gas turbine 


It should be made clear at the start 
that the term gas turbine is generic, 
and refers to the internal combustion 
gas turbine. 


As such, the gas turbine came into 
the public eye during the latter part 
of World War II as a power unit for 
high speed combat airplanes. In the 
early days the main difficulties were 
obtaining a compressor of  suffi- 
ciently high efficiency that would 
not consume too large a proportion 
of the power generated by the gas 
turbine, and obtaining materials that 
could withstand the high tempera- 
ture of combustion (1300 to 1400°F). 


Since World War II the gas tur- 
bine has undergone considerable de- 
velopment and reached a compara- 
tively high standard of performance 
and reliability in the aviation field. 
In the industrial field, however, its 
development and utilization has not 
been quite so rapid. 


The reason for this is apparent if 
we review the priorities of the main 
factors affecting the gas turbine de- 
sign. 

(See top of page 61, column 2 and $3) 


The Canadian Journal of Chemical Engineering, April, 1958 


Se ee 


a oe 2am 2 = Ge Gee os 


a-_m -= 2 es em 


ee ee 


TT Oem eR mM emeO Fn OUlklOF 


oonooofoe!-_ 


as 
Or 


bt 
ill 


8 


a the aviation field the gas turbine 
has offered a performance so much 
in advance of other engines of the 
same weight and volume that its 
otentialities could not be ignored. 
It thus fulfilled the main combat re- 
quirements and automatically ob- 
tained the capital necessary for 
research and development. Conse- 
quently, the aircraft gas turbine 
reached a high standard of perform- 
ance and of utilization. (See com- 
parative weights of gas turbines and 
piston engines in Figure 1). 

In the industrial field the ability 
to produce a high power output w ith 
a low weight and space requirement 
has a limited application other than 
for automobiles, speed boats and rail- 
way traction. Accordingly, the selec- 
tion of a gas turbine driver for large 
process compressors and refrigera- 
tion compressors becomes mostly a 
matter of economics. As the first cost 
of the gas turbine is high and its 
efficiency is low, it is generally a 
matter of arithmetic to determine 
whether or not a particular installa- 
tion will be economical, which de- 
pends mostly on fuel cost. 

There are four main types of gas 
turbine: 

(a) Simple cycle, single shaft. 

(b) Simple cycle, double shaft. 

(c) Simple cycle, with regenerator. 

(d) Compound cycle with  inter- 
cooler and regenerator. 


In the simple cycle, single shaft 
machine shown in Figure 2, fresh 
air is drawn into the compressor and 
compressed and discharged into the 
gas turbine. Fuel is pumped into the 
air after compression in the compres- 
sor and prior to entering the gas 
turbine. In the gas turbine, combus- 
tion takes place and the hot gases 
expand through the turbine blading 
before passing out to the atmosphere. 
In the single shaft machine the gas 
turbine is directly connected to the 
compressor shaft. As the compressor 
gives high efficiency only at full 
speed, the single shaft machine effi- 
ciency is very low at part loads be- 
cause at low speeds the internal work 
required to drive the compressor 
represents the major portion of the 
power generated by the gas turbine. 

In the two shaft or “free turbine” 
type of machine there are two tur- 
bines both gas driven. The first is 
directly connected to the compressor 
which it drives, the second turbine 
provides the power for the load. 
Thus, it is possible to run the com- 
pressor at optimum efficiency speed 
while the speed of the power turbine 
can be varied to suit load require- 
ments. An additional advantage of 


Gas turbine design factors in order of priorities 


Aircraft Gas Turbine 


1. Performance and reliability. 

2. Weight and space. 

3. Operating economy — fuel and 
maintenance economy. 

Initial cost. 


Industrial Gas Turbine 
Operating costs, fuel equipment. 
Initial capital cost. 

Reliability and performance. 


- wen — 


Weight and volume. 4. 


WEIGHT COMPARISONS 
GAS TURBINE 


AIRCRAFT 
GASOLINE 


AUTOMOTIVE 
GASOLINE 


MARINE DIESEL 





VOLUME COMPARISONS 


GAS TURBINE 





AIRCRAFT 17.6 CU. FT. 


GASOLINE 


AUTOMOTIVE 
GASOLINE 


MARINE DIESEL 


Figure 1—Weight and volume comparisons of gas turbine and piston engines 
of equal rating. 


4 IN THE COMBUSTION CHAMBER THE COMPRESSED 
AIR 1S HEATED BY THE BURNING FUEL 


5 HoT GASES ARE 
PASSED THROUGH 
THE TURBINE BLADES, 
THEY CAUSE 

THE TURBINE 

TO ROTATE 


FUEL IS PUMPED 
INTO THE COMBUSTION 
CHAMBER 





POWER 
TAKE OFF 
COUPLING 


THE AIR IS 

COMPRESSED BY THE 

COMPRESSOR, AND IT 

FLOWS TO THE 

COMBUSTION CHAMBER 6. EXHAUST GASES ARE 
PASSED OUT INTO THE 


ATMOSPHERE 


| FRESH AIR 18 DRAWN 
* INTO THE COMPRESSOR 
INTAKE 
Figure 2—Schematic representation of a simple open-cycle gas turbine for 
mechanical shaft drive. 
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PERCENT OF RATED LOAD 


Figure 3—Approximate fuel rate curves for steam turbines, diesels, and simple 


gas turbines (at various load levels). 


TABLE 


1 


lypIcAL EFFICIENCIES AND HORSEPOWERS OF VARIOUS THERMAL DRIVERS” 


lype of Power Plant 


Diesel Engine 

Large Steam ‘Turbine (Condensing ) 
Medium Steam Turbine 

Small Steam Turbine. . . 

Steam Locomotive 

Gas Engine (Reciprocating ) 
\utomotive Gasoline Engines 

Best Gas Turbine (Complex) 
Simple Open-Cycle Gas Turbine 


Small Simple Open-Cycle Gas Turbine. 


this arrangement is that the power 
turbine has a high starting torque. 
This starting torque is at a maximum 
at a stand-still. This feature makes 
the two shaft turbine ideally suited 
for starting heavy loads. 

The third type of gas turbine is 
the simple cycle with regenerator. 
\ simple cy cle may employ one or 
two turbines as in (a) or (b), but 
the regenerator, or recuperator, takes 
the form of a heat exchanger which 
heats the discharge air from the 
compressor to the combustor by 
utilizing the heat of the exhaust gases 
from the gas turbine. By this means, 
it is possible to increase the overall 
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Typical Size Ranges 
Thermal | in 
Efficiency Horsepower 


30-34% 20-7000 


28-35% | Up to 350,000 
20-25% | Up to 27,000 
16-20% | Upto 1,000 
6-10% | Up to 2,000 

25% | Upto 2,500 

22% | 40-300 
30.5% | 36,000 
14-18% 1,000-25 ,000 
6-12% 50-500 


efficiency from 16% to 23%. Typi- 
cal efficiencies for various types of 
gas turbine and for other ‘internal 
combustion engines are given in 
Table 1. Approximate fuel rates are 
shown in Figure 3. 

For large efficient industrial in- 
stallations the compound cycle with 
intercooler and regenerator is em- 
ployed. Here air compression is car- 
ried out in two separate compressors 
in two stages with an intermediate 
cooler in between, thus improving 
the compression and overall dynamic 
efficiency. The efficiency is further 
improved by heating the compressed 
air by means of the exhaust gases 


from the gas turbine. In some cases, 
as for example with a peak load tur- 
bine, the heat exchanger or recupera- 
tor is omitted. Two stage compres- 
sion with an intermediate cooler not 
only increases the efficiency of the 
unit but also gives a_ substantial 
power increase for the same volume 
of suction air. This reduces the cost 
per installed kilowatt. 

Here are some typical figures re- 
lating to a two shaft, non-recupera- 
tive gas turbine. The particular unit 
considered is a Ruston Hornsby type 
TA gas turbine. This unit is rated 
at 1,260 continuous horsepower with 
an efficiency of 15.9% and a maxi- 
mum overload rating of 1,410 horse- 
power at 16.2% efficiency. The 
normal speed of the power turbine is 
6,000 RPM. Inlet air conditions are 
assumed to be 14.7 p.s.i.a. 65°F. After 
compression in the compressor these 
conditions become 61 p.s.i.a. 380°F. 
After the compressor the compressed 
air goes to the combustion chamber 
at a pressure of approximately 60 
p.s.i.a. and a temperature of 1,340°F. 
is achieved due to the combustion of 
the fuel. After passing through the 
first turbine, that is the turbine driv- 
ing the compressor, the gas condi- 
tions are 24 p.s.i.a. and 1,020°F. at 
which conditions they enter the 
power turbine wheel to emerge from 
there at 14.7 p.s.i.a. at 845°F. 

The efficiency of the gas turbine, 
16.2% without regenerator or 23.1% 
with regenerator under maximum 
rating or overload conditions, is low 
compared with 22% for the auto- 
motive gasoline engine, 25% for the 
reciprocating gas engine, 20-25% for 
a small steam turbine, 28-35% for 
large steam condensing turbine and 
30-34% for a diesel engine. How- 
ever, the heat immediately lost to 
shaft horsepower is recoverable. 

There are three alternative forms 
of heat recovery. The first is to use 
the heat of the turbine exhaust gases 
in a heat exchanger or regenerator 
to heat the compressed air. The 
second is to use all or part of the 
exhaust gas heat in waste heat boilers. 
These waste heat boilers may be en- 
tirely fired from the exhaust gas 
from the turbine or may have sup- 
plementary fuel firing as_ well. 
Thirdly, the exhaust gases may be 
used to heat air for industrial pur- 
poses. 

Whereas, in England, Ruston and 
Hornsby have standardized on a 
unit of some 1,400 horsepower; in 
the States, the General Electric Com- 
pany has produced a unit of 5,000 
horsepower. This unit is a simple 
cycle single shaft gas turbine prim- 
arily intended for generation or 
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mechanical drive at fixed speed and 
is now being manufactured in the 
U.S.A. General Electric also make a 
two shaft regenerative cycle gas tur- 
bine. 

The Westinghouse Company have 
standardized on an 1,800 horsepower 
gas turbine plant of the simple open 
cvcle type; they also produce un‘ts 
of 5,000 horsepower and 7,000 horse- 
power. The major component parts, 
the axial flow air compressor, the 
combustion system, and gas turbine 
are arranged in line and mounted on 
a fabricated steel bed-plate. The 
complete unit including all acces- 
sories is factory assembled and can 
be shipped in one package. It oper- 
ates at a speed of 8,750 RPM. The 
temperature of the hot gases after 
the fuel has been admitted is 1,250°F. 
About two-thirds of the power de- 
veloped by the gas turbine is used 
to drive the axial flow air compressor. 
The remainder is available as useful 
power. This unit gives approximately 
119,000 Ibs./hr. of exhaust gas at a 
temperature of 780°F. for waste heat 
recovery. 

Theoretically, the gas turbine can 
be fired with almost any fuel. No 
particular octane nor cetane number 
is required in order to obtain satis- 
factory performance. The preferred 
fuel is natural gas although sewage 
gas, blast furnace gas and similar 
gases have been used successfu'ly. 
Generally speaking, any clean gas 
fuel may be used. Alternative'y, 
gasoline, kerosene, gas oil or diesel 
oil and distilled fuel oils can be 
burned without trouble in a gas tur- 
bine. 


Residual fuel oils have been used 
but their natural impurities often 
cause serious problems. Impurities 
such as vanadium and sodium, some- 
times in the form of the pentoxide in 
the case of vanadium, are frequent!y 
encountered in South American 
crude oils. These impurities melt at 
about 1,240°F. forming slags which, 
in the case of vanadium, have an ex- 
tremely corrosive effect upon the 
high temperature alloys used for the 
turbine blading and the combustion 
chamber parts. As many modern gas 
turbines are commonly designed to 
operate with inlet turbine tempera- 
tures in the 1,400-1,500°F. range, the 
corrosive content of residual fuels 
can cause extensive damage to the 
turbine. In view of the fact that 
residual fuels are often available at 
half the cost of distillate oils, inten- 
sive research is being carried out in 
Switzerland, Britain, and the United 
States to develop suitable methods 
of burning such residual oils. In 


TABLE 2 


Gas TURBINE 
| 


Simple-Cycle 


CHARACTERISTICS Data‘ 


Simple-Cycle 


Regenerative 








| 
Single-Shaft Two-Shaft | Cycle 
Rating (80°F and 1000 ft. 
PAPRTRIIG: NEMS es cos css See ee 3500 kw 3500 kw 500 kw 
*Fuel Rate (BTU/KWH)......... 23,900 22,900 16,000 

Over-all Eiicieney. oc. sic .cce cee 14.3% 14.9% 21.4% 
Maximum Overload Rating (Low 

Ambient Temperature)...... 5000 kw 5000 kw 5000 kw 
Cooling Water Required (GPM)...} 300 400 400 
Lubrication oil tank capacity 

CORIO «i.e tie ncn fowen wee aaee 1100 900 900 
Dimensions 

PRE 9.35, tiene Katee 49 ft. 5S ft. 53 ft. 

WENN oso oe) oe eh a wees 10 ft. 11 ft. 14 ft. 

RO Pi ano accor eat eee 10 ft. 10 ft. 10 ft. 
Weight (Ibs.).............e.00e- 107,000 130,000 380,000 

**Maximum Weight after 

Installation (]b8.).......0%6o6%%. | 30,000 30,000 30,000 

COR ies a ede oe rer en ren Natural gas or Natural gas or | Natural gas or 
Distillate oil | Distillate oil Distillate oil 

Air Requirements (80°F. and 

1000 ft. Altitude Max.)........ 98 lbs. /sec 98 Ibs./sec 98 Ibs./sec. 
Standard Generator..............| 3500 kw .8 pf | 3500 kw .8 pf | 3500 kw .8 pf 

| 


Oversize Generator.......<....5.- 


***Top Mounted Air Cooler (water 





4000 kw .8 pf 


4000 kw .8 pf 





4000 kw .8 pf 


temperature standard)......... | 85°F. | 85°F. 85°F. 
weeP Standard VORAGES: ss. c ccc a ws | 2400/4160 | 2400/4160 2400/4160 

Direct-connected Exciter......... | 25 kw | 25 kw 25 kw 
Slip Ring Starting Motors........ 250 hp 250 hp 250 hp 
Gas Turbine Starting and s ; os 

Supervisory Control........... included included | included 
Turbine Auxiliary Control Panel. .| included included | included 

| 


All units can be rated 5000 hp for Mechanical-drive applications. 


*Based on higher heating value of Natural Gas fuel of 1,100 BTU 


value 1000 BTU/cu. ft.). 


cu. ft. (Lower heating 


**Based on skidding the large stationary parts in place such as generator armature and 


regenerators. 


*** Air coolers for 90 and 95°F. temperature can be furnished. 
****Other voltages such as 6,900, 11,500 and 13,800 volts can be furnished. 


parallel with this research, attempts 
are being made to develop solid fuel 
fired gas turbine utilizing coal, peat 
and sawdust. 

While discussing fuel, it shou'd be 
mentioned that the gas turbine is 
very sensitive to temperature of the 
ambient air; a 40°F. decrease in inlet 
air temperature will result in a 20% 
gain in output power. 

The gas turbine like the centri- 
fugal compressor has few rotating 
parts; consequently, mechanical wear 
is reduced to a minimum and vibra- 
tion is less than in the case of re- 
ciprocating engine. On account of 
the high operating temperatures, 
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however, the blading of the gas tur- 
bine is subject to considerable cor- 
especially when using fuels 
wh-ch do not burn cleanly. 

Another feature of the gas turbine 


rosion, 


is the 


very 


small consumption of 


lubricating oil, especially when com- 
pared with a reciprocating engine. 
With regard to maintenance costs 


the General Electric Company have 
given a figure of 0.003 cents/kwhr. 
of generated electricity for a 4,000 
kw machine running for a period of 
3,500 hours at the Belle Isle station 
of the Oklahoma Gas and Electric 
Company. A comparable figure for 
steam plant maintenance is 0.04 
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TABLE 3 


\ CoMPARISON BETWEEN ELECTRIC Motor AND STE AM TURBINE Driv ES 


Electric Back Pressure Steam Turbine 

Motor 

(Base 600 psig 900 psig 1300 psig 

Case) inlet inlet inlet 
Steam rate- he. /hp hr. 26.2 20.5 18.4 
Investment—$/hp 118 157 206 262 
Net over base case—$/hp 39 88 | 144 
Annual charges—$/hp yr. ed 23 | 27.6 32.9 
Gross Savings over 

base case—$/hp yr. 49.7 45.1 39.8 

% Net return 57 23 12 


cents/kwhr. while for modern diesel 
engines it would be 0.2 and 0.3 
cents/kwhr. of generated electricity. 

A. further advantage of the gas 
turbine is that cooling water require- 
ments can be reduced to zero. This 
is particularly advantageous in loca- 
tions where gas turbines are used in 
the desert on pipe line installations 
where water is not obtainable. It is 
of particular interest to oil com- 
panies where present day thinking in 
favorable localities is in terms of air- 
cooling rather than of water cooling. 

Compared with the steam conden- 
sing turbine, the gas turbine requires 
smaller foundations. 

In two of the cases studied, we 
have considered waste heat boilers 
for process steam generators. About 
6-7 pounds of steam per horsepower 
hour of 150 p.s.i.g. steam can be gen- 
erated for about a 5-10% increase in 
the gas turbine investment cost. 

One gas turbine application of in- 
terest is the utilization of a gas tur- 
bine to provide combustion air for 
fluid cat. cracker regenerators. With 
the open cycle type of gas turbine, 
it is customary to supply some six 
or seven times the quantity of air 
required d for combustion in order to 
maintain a gas temperature not ex- 
cessively hot for the turbine blades. 
The idea was conceived by a former 
member* of our utilities group to 
utilize the exhaust gas, including un- 
burned air from a gas turbine, to 
provide combustion air for the re- 
generator in a fluid cat. cracker. 
Originally, it was proposed that the 
gas turbine be so designed that it 
would not have any excess power 
available at the shaft, but further 
developments allowed for excess 
power being available to drive the 
cat. cracker gas compressor. Some 
characteristics of the three types of 


ey. Guthrie. 
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gas turbines are shown in Table 2 
For stationary industrial installa- 
tions, it would appear that at present 
the gas turbine can only be com- 
petitive where cheap or non-market- 
able fuel is readily available in the 
form suited to the gas turbine. As 
special metals for gas turbine blade 
construction become available and 
design techniques are developed to 
overcome high temperature difficul- 
ties in the blading and combustion 
chamber, the gas turbine will un- 
doubtedly be used increasingly. 


3. Steam turbines 

Most pumps and compressors in a 
refinery are driven by electric 
motors, directly connected whenever 
possible. In other cases, steam is used 
on account of its reliability and 
flexibility; also, to provide exhaust 
steam for tank heating. 

Steam turbines as used in a refin- 
ery can be divided into two broad 
groups. 

(i) Single stage mechanical drive 
steam turbines up to 1,000 horse- 
power. 

(ii) Multi stage units of 1,000 horse- 
power or more for the larger 
and more important drivers such 
as cat. cracker blowers and com- 
pressors and the process and re- 
frigerator gas compressors for 
the Petrochemical Plant. The 
power required for these tur- 
bines varies from 1,000 horse- 
power to 8,000 horsepower with 
speeds from 3,000 RPM to 9,000 
RPM. 

For such duties, the steam turbine 
unlike the induction or synchronous 
motor provides variable speed offer- 
ing greater operational flexibility and 
the ability to operate under changed 


design conditions such as different 
molecular weight in the case of gas 


compressors; also, to operate at part 
load without recycling the fluid be- 
ing compressed. 

Reliability in general is very good, 
availability up to 98% and even 99% 
being quoted by some users. Gen- 


erally, back pressure turbines are 
used, which eliminates any opera- 
tional difficulties with condensers 


and may provide cheap power ac- 
cording to the amount of exhaust 
steam that can be used. 

In cases where back pressure tur- 
bines are installed, desuperheaters 
operate in parallel to maintain the 
supply of yard steam when the back 
pressure turbines are not running, 
The boiler investment cost, there- 
fore, just covers the incremental cost 
of the high pressure boilers above 
that of the low pressure boilers. 

Typical costs for back pressure 
turbines exhausting to 110 p-s.i.g. 
compared with electric motors of 
92% efficiency operating on pur- 
chased power are given below in 
Table 3. 


In Table 3, boiler investments are 
based on a capacity of 200,000 lbs./ 
hr. with a boiler etficiency of 86.5%. 
Figures tabulated show the difference 
between the cost of a 150 p.s.ig. 
boiler plus the steam distribution sys- 
tem and the cost of the high pressure 
boiler plus its steam distribution sys- 
tem for the steam rates quoted. 
Water treatment costs are order of 
magnitude investments plus operat- 
ing chemical costs which can vary 
considerably. 


Operating labor is incremented 
over the base case. Forty per cent 
excess boiler capacity is included in 
the investment costs to allow for 
maximum 110 p.s.i.g. steam require- 
ments plus some “spare capacity”, 
although at present we do not pro- 
vide for spare boiler capacity. The 
turbine horsepowers are, of course, 
limited by the minimum 110 p-S.i.g. 
requirements in summer. 

The horsepower considered was 
1,800 at 100% capacity and 100% 
serviceability. Depreciation was 5%, 
and income tax 55% of gross savings. 


The above figures (courtesy Esso 
Engineering) confirm that 600 p.s.i.g. 
is at present the highest economical 
operating pressure with the condi- 
tions given. 

Recent experience with 600 p.s.i.g. 
operation (and to some extent of 400 
p-s.i.g. also) has shown that the foul- 
ing of turbine blades can be experi- 
enced at this pressure due to silica 
deposits. 

At Sarnia, a raw water content of 
7 PPM silica. (as SiO,), is reduced 
to 2 PPM after treatment in a hot 
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lime soda softener, resulting in a 
figure of 18 PPM in the boiler 
druins at a concentration of 9 times. 
A'though steam purity has _ been 
maintained at less than 1 PPM, loss 
of power on the cat. cracker blower 
turbine has resulted after periods of 
6 to 8 weeks. As a result, turbine 
washing facilities have been installed. 
Subsequent experience has shown, 
however, that the turbine blade foul- 
ing may not have been entirely due 
to silica, but to another impurity 
which helped the silica to adhere to 
the turbine blades. 

Although “offload” washing equip- 
ment is cheaper to install and safer 
to operate than “onload” equipment, 
the latter is preferred to avoid down- 
time of the unit. Experience has re- 
sulted in a policy of sparing turbine 
rotors in the case of “on load” 
stallations due to the hazards of ex- 
cess thrust and damage to the turbine 
thrust gear. This in turn entails a 
considerable increase in investment 


cost. (See Table 4). 


Recent experience has also shown 


COMPRESSOR DRIVERS FOR PETROCHEMICAL PROJECT INVESTMENTS 


Case: I 


Total cost of boilers, auxiliaries 
and utilities plant. . 


Distribution system at site 57,500 | 57,500 86,200 57,500 | 
Dest iperheating station 
(150,000 Ibs./hr.). . 50,000 50,000 
12” dia. 4000’ hp steam line 
from Plant No. 2. 120,000 
Total cost of waste heat boilers 
and auxiliaries—installed...... 267,000 
Electric Power—Main Sub. and 
distribution 769,000 709,000 191,000 191,000 
Electric Motors... . 232,000 232,000 
Steam Turbine. 80,000 80,000 
Steam Turbine Washing 
Facilities. . 75,000 75,000 | 
Gas Turbines: (a) 8500 hp 568,500 583,500 | 
(b) 6000 hp - 408,000 423,000 
(c) 1800 hp. ~ 216,000 | 231,000 | 
Installation of Drivers... - 85,000 85,000 | 
| 
Gas Turbines Fuel System. 25,000 25,000 | 
Spare Rotors. 30,000 30,000 360,000 360,000 
Total.... $2,491,000 | $2,358, 000° $2,856,700 | $2, 931, 000 


*Based on $6.5 per pound per hour. 
tBased on $7.5 per pound per hour. 


$*1,137,500 


that the number of auxiliaries and 
safety devices apparently needed on 
large back pressure turbines can re- 
duce reliability appreciably. It would 
appear w iser to simplify the installa- 
tion as much as possible and to de- 
sign for rugged construction rather 
than high efficiency, the only merit 
of which is to reduce the investment 
cost of the high pressure boilers. 
An efficient turbine requires several 
stages and hence a long turbine shaft 
which is subject to thermal shock. 
This resulted in “bowing” and out 
of balance due to a sudden shutdown 
in the case of a 600 p.s.i.g. turbine 
driving a main air blower at Halifax 


Refinery. 


In the larger Imperial Oil refin- 
eries, steam is used at three different 
pressures: 


600 p.s.i.g. 
125 p.s.i.g. 
15 p.s.i.g. 


High Pressure 
Yard Steam 
Exhaust Steam 


High pressure steam is used for a 
few large turbines which exhaust 
into the yard steam system. 


TABLE 4 


| 


II IV V 


$71,125,000 | $ *650,000 | $ *975,000 
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Yard steam supplies most of the 
process requirements, such as pump 
turbines, steam for vacuum jets, par- 
tial pressures, stripping steam, purg- 
ing, aeration, heating in reactors and 
for jackets, and for line blowing 
aut start-up. 

Yard steam is also used for tank 
heating where a temperature over 
150°F. is required. 

Exhaust steam is used for tank 
and building heating, for snuffing 
steam and for gas freeing tanks. In 
the utilities plant, clean 15 psig. 
steam from turbine exhausts is used 
for heating boiler feed water in the 
deaerator. 


We thus have a cascade system 
where steam is generated first at 600 
p.s.i.g. to supply the one or two large 
back pressure turbines driving criti- 
cal compressors and blowers. The 
exhaust steam from these large tur- 
bines, which is clean (i. free of 
lubricating oil), discharges to the 
yard system. At this level, some of 
the steam is recoverable, such as the 


VI 
975,000 
57,500 
Only in Case I will there be 
a steam tie-in with Sarnia 
Refinery 
797,000 | Estimate based on $125 per 
KVA 
| 
296,000 | Motors are weather 
protected type 
- Motor and steam turbine 
assumed the same. Fig. 
| shown is incremental only 
for gas turbines 
$2,125,500 
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exhaust steam from 125 p.s.i.g. tur- 
bines, and some, such as stripping 
steam, is not. In parallel with the 600 
psig. to 125 p.s.i.g. steam flow, 
there are desuperheaters which pro- 
vide continuity of 125 p.S.L.g. steam 
supply should the high pressure tur- 
bines fail, and additi ‘onal 125 p-S.L.g. 
steam if required. Usually the total 
desuperheating capacity numerically 
equals the high pressure boiler 
capacity. 

The exhaust system has “make-up” 
connections at suitable points con- 
sisting of pressure reducing valves 
supp'ied with 125 p.s.i.g. steam. 

With such a system, the horse- 
power that can be provi ided by high 
pressure turbines is limited to the 
corresponding amount of exhaust 
steam from these turbines which can 
be used at the 125 p.s.ig. level — a 
factor in the selection of drivers. 
Alo, while the supply and demand 
of steam can be balanced at the high 
pressure level (by control of the 
hi gh pressure boilers), and at the 125 
p.s.i.g. level by means of the de- 
superheaters, it is not always possible 
to absorb all the exhaust steam, par- 
ticularly in the summer when build- 
ing and tank heating requirements 
are small. 

Any prolonged excess of exhaust 
steam, apart from being wasteful, 
completely upsets the economics of 
" back pressure turbines. 


Figures for the total heat of the 


steam above 32°F. as taken from the 
team tables are: 


When considering the overall 
economics of the system, the total 
heat per pound is corrected to allow 
for the fact that steam is raised from 
water at 50°F. and not from water 
at 32°F., where 50°F. is the average 
temperature of the raw water 
throughout the year at the refinery 


High Pressure 
600 p.s.i.g. and 750°F. 
Yard steam 
125 p.s.i.g. and 400°F.* 1221.8 _ 
Exhaust steam 
15 p.s.i.g. and 300°F.* 1189.2 - 
*(slightly superheated due to high pressure 
turbines and desuperheater). 


1378.2 BTU/Ib. 





TABLE 5 


COMPRESSOR DRIVERS FOR PETROCHEMICAL PRoJECT CASES STUDIED 
Case I 
(a) One 175,000 lbs./hour boiler located at Plant No. 2*; steam at 650 p.s.i.g. and 750°F. 
(b) One 3,200 horsepower back pressure steam turbine driving propylene compressor. 
(c) Total 11,600 horsepower of electric motors for all other compressors. 
(d) Four thousand (4,000) ft. of 12” dia. high pressure steam line from Plant No. 2 to 
petrochemical area. 


Case II 

(a) Two 75,000 lbs./hour boilers located at petrochemical area; 
and 750°F 

(b) One 3,200 horsepower steam turbine driving propylene compressor. 

(c) Total 11,600 horsepower of electric motors for all other compressors. 


steam at 625 p.s.i.g. 


Case IV 

(a) Two 50,000 Ibs./hour boilers located at petrochemical area; steam at 150 p.s.i.g. 
dry sat. 

(b) Three waste heat boilers generating a total of 89,000 Ibs./hour steam at 150 p-s.i.g., 
dry sat. 


(c) Three gas turbines (8500, 6000 and 1800 horsepower) driving all compressors. 


Case V 
(a) Two 75,000 Ibs./hour boilers located at petrochemical area; steam at 150 p.s.i.g., 
dry sat. 


(b) Three gas turbines (8500, 6000 and 1800 horsepower) driving all compressors. 
No waste heat boilers. 


* Plant No. 2 is some 4,000 feet west of the petrochemical main site. 
Case VI 
(a) Two 75,000 Ibs./hour boilers located at petrochemical area; steam at 150 p.s.i.g., 


dry sat. 


(b) All compressors driven by electric motors—total 14,800 horsepower. 


TABLE 6 


CoMPRESSOR DRIVERS FOR PETROCHEMICAL PROJECT 


Case: 


Tot: il larvieeniaelt Costs. 


Fuel Fired in Boilers (Bunker Oil)*.... 


Fuel Fired in Gas Turbines (Diesel Oil) 


Fuel Fired in Gas Turbines (Gas and Bunker 


Price). . 
Power for Drivers only. . 
Operators (3 shifts per day)... 
Maintenance (4%)... 
Depreciation (5°) 


Total Manufacturing Costs. . 


Sitersnce i in levestunent over Case VI. 


Savings in Manufacturing cost over Case VI.... 


Income Tax (47%) 


Net Savings (53%) ..........5. 


% Net Return on Difference in _Investment 


OPERATING Cost AND RETURNS ON INVESTMENT 





I II IV Vv VI 

$2,491,000 | $2,358,000 $2,856,700 | $2,931,000 | $2,125,500 
}. (a) |. (b)__i 

935,000 | 935,000 | $ 294,000 | $ 294,000} 798,000} 798,000 

$1,478, C00 — - 

atid - | - | - | 868, 000 | - - 
640, 000 640,000 | ~ * | - | 817,000 
25,000 75,000 100,000 100,000 | 75,000} 75,000 
100 ,000 94,000 114,000 114,000} 118,000} 85,000 
125,000 118,000 143,000 | 143,000} 147,000 | __ 106, 000 


$1,881,000 


$1,825,000 | $1,862,000 | $2,129,000 | $1,519,000 $2,616,000 | 


$ 365,500 | $ 233,000. ) $ 


731,200 

56,000 19,000 362,000 

iis 26, 300 8,940 170,000 
29,700 10,080 192 ,000 | 
er 8.1% 4.3% ! 26% | | 











* or Refinery Gas if available at Bunker prices. 
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TABLE 7 


COMPRESSOR DRIVERS FOR PETROCHEMICAL PROJECT 
ESTIMATED STEAM, POWER AND WATER REQUIREMENTS 


| 125 psig 


Unit Steam | 
| (Ibs./hr.) 
Ste am Cracker +Lighe Ends _ | 78,000 _ 
Ethylene Recovery | - 
Clay Plant | 2,800 
Butadiene Extraction 4,000 
Heating Oil Hydrofiner | 250 
Total Process | 85, 050 | 
Buildings | 20,000 
Tank Heating | 10,000 
Utilities (18.5%) | 26,150 
Offsite | + 
Total “141,200 | 





Water USGPM 


Drive 
H orse Power 
Power (KW) Cooling Fresh 
t =30°F. Process 
7,500 | 820 14,800 100 
1,500 | 150 2,800 | 20 
3,200 | 35 
~ 1. - 2 2 
| 
600 | 210 400 | 50 
= | is | 210 
14,800 | 1,230 19,260 | 170 
| 
| | | 
| 
| 
| 1,700 
| 2,930 


t 30 lbs./Horse Power/Hour and assuming 141,200 lbs./Hour maximum Horse Power 


available from back pressure turbine is 141,200 


30 


= 4,700. 


At 30 lbs./Horse Power/Hour and assuming 85.050 + 18.5% max. Horse Power available 


from back pressure turbine is 101,000 
30 


in question. The poy figures then 
are reduced by 18° and become: 


for high pressure steam 1360.2 BTU's/Ib. 


‘* yard steam 1203.8 Pe 
“* exhaust steam 1171.2 e 
(When calculating efficiencies we 


have to distinguish between “boiler 
efficiency” when the inlet feed water 
temperature at the boiler would be 
used (usually 220°F.), and “overall 
or system efficiency” when the yard 
raw water temperature (50°F.) 
would be used). 

In the attached tables, fuel is 
priced at so much per barrel of 6.485 
million BTU’s. Provided there is no 
surplus and therefore wastage of 
steam, it would seem reasonable to 
charge high pressure turbines with 
1360.2 - 1203.8 

1360.2 
of raising high pressure steam, and 
1203.8 - 1171.2 

1203.8 


of the cost per pound 


in the case of low 


pressure turbines. This procedure 
will generally show economics in 
favor of back pressure turbines and 
is, of course, the reason for installing 
“topping” or back pressure turbines. 
The economics will be completely 
upset, however, if the remaining heat 
in the exhaust steam is wasted as it 
represents so large a proportion of 
the total heat in the high pressure 
steam. 


= 3,300. 


4. Electric motors 


In 1954, the petroleum industry in 
the U.S.A. consumed over twelve 
billion K.W. hours annually of elec- 
trical energy. Refining accounted for 
eight billion K.W. hours of this. In 
1940, the total horsepower hours 
used per barrel in refining was 
divided 44% amongst electric drives 
and 56% amongst steam drives. 

In 1954, the figures were: 


For small For large 


plants plants 
Electric 100% 40% 
Steam 0% 60% 


These statistics show the present day 
trend. The figures for large plants 
are not representative as they apply 
in most cases to one or two large 
units per refinery (motors over 1,000 
horsepower). 

Electric drive is preferred for low 
investment cost, low maintenance 
and simplicity in operation. Its pos- 
sible disadvantages are the depend- 
ability of the supply source, and the 
particular installation. 

The Squirrel cage induction motor 
is a mechanically robust and reliable 
machine, while the synchronous 
motor can offer improved economics 
by virtue of power factor improve- 
ment and high efficiency. Both are 
no more reliable than the source of 
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power they consume, and both lack 
speed variation. The induction motor 
has the disadv: antage of the amount 
of current and time required to 
reach full speed which may be limit- 
ing factors where the supply source 
has insufficient back- up capacity. Ac- 
celeration times of 20 to 30 seconds 
are not uncommon 

the moment of inertia 


resulting from 


(WR2) of the 


motor (2,500 Ib. ft.2) and of the 
centrifugal compressor (13,000 Ib. 
ft.2). 


Although the synchronous motor 
can be used in cases of high inertia 
drives, it is difficult to ‘apply in 
hazardous atmospheres. 

Imperial Oil has two motor driven 
cat. cracker blowers at loco and 
Winnipeg and a gas compressor unit 
at Montreal. The availability and re- 
liability of these installations has 
compared very favorably with other 
steam driven units. A table of known 
electrical drives for cat. cracker main 
air blowers is given below: 

Atlantic Refining Company (2) 

Philadelphia 
Bay Petroleum Company 

Denver, Colorado 
Ben Franklin Petroleum Company 

Ardmore, Oklahoma 
British American Oil Company 

Toronto, Ontario 
Eldorado Refinery Company 

Lamont, Illinois 
Imperial Oil Limited 

Winnipeg, Manitoba 
Imperial Oil Limited 

loco, British Columbia 
Kanotex Refining Company 

Arkansas City 
Northwestern Refining Company 

. Paul Park, Minn. 
Rock Island Refinery 

Indianapolis ; 
Standard Oil Company 

Luna, Ohio 
Standard Vacuum Refining Company 

Durban, South Africa 
Canadian Petrofina Limited 

Montreal, Quebec 


Two distinct advantages of the 


electric motor are its operational 
simplicity it is usually push-button 
controlled — and the fact that it does 


not require an operator in attendance. 
For the petrochemical project, the 
choice of electrical motor drives de- 
pends almost entirely on economics, 
as power rates and continuity of sup- 
ply are favorable, and ample capa- 
city can be made available in the 
Sarnia area for the horsepowers we 
require. 
5. Comparison of the various cases 
considered* 


During preliminary studies, a num- 


®*See Tables 4-8. 
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ber of cases were considered based 
on different locations for the utilities 
plant and different types of drivers. 

It appears that for a certain size 
refinery — perhaps 75,000 barrels per 
day — it is economical to centralize 
the utilities plants, but beyond that 
size, it is better to decentralize. Ac- 
cordingly, in some of the cases 
studied the petrochemical utilities 
centre is located at the existing 
Sarnia Refinery No. 2 Plant; in 
others, it is located at the new petro- 
chemical plant site. The most attrac- 
tive of the several cases studied are 
described on Table 5; Table 4 gives 
investment costs. 

Table 6 gives operating costs and 
returns on investments based on the 
estimated utility costs and fuel costs 
on Tables 7 and 8. 

The estimate includes feeders from 
the Hydro substation, the main sub- 
station and feeders to the compressor 
motors. It does not include unit sub- 
stations nor feeders to other large 
motors. Electric power is available 
at 27 KV and the supply is con- 
sidered as reliable as grid steam 
plants. No difficulties were expected 
as regards starting large electric 
motors. 


6. Conclusions 

(i) Cases I and II show a small re- 
turn on the incremental invest- 
ment, but the return is not large 
enough to warrant the increase 
of investment. The small return 
is due to the fact that the horse- 
power of the back pressure tur- 
bines is limited by the 125 
p-s.i.g. steam requirements to 
3,200 horsepower. 

(ii) If refinery gas can be made 
available, and at a price of $2.80 
per fuel oil equivalent barrel, 
case IVb shows a return of 26% 
on the incremental investment. 
This is the best economic case. 
Before recommending gas tur- 
bines to management, we must 
check: 

(a) Availability of gas fuel. 

(b) The present inability of gas 
turbine manufacturers to 
guarantee successful opera- 
tion on bunker and other 
residual oils, as alternative 
fuels. 

(c) The uncertain, and in some 
cases, long delivery time.* 

(iii) If fuel gas is not available in 
sufficient quantity for gas tur- 
bines, Case VI, the “All Elec- 
tric” Case, offers the lowest 
capital investment, and the ad- 
vantage that the cost of electric 
power is reasonably steady and 





®See Table 9. 
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TABLE 8 


COMPRESSOR DRIVERS FOR PETROCHEMICAL PROJECT 
_ Cost AND > BTU VALUE oF Fu ELS 





Dollars BT U's BT ‘U" s BTU's Dollars 
Fuel we - — nae . 
Barrel end cu u. ft - Dollar A 000, 000 BTU 
Bunker Oil *2.80 6,275,000 240, 000° 0.447 4 
Diesel Oil *4.50 5,930,000 1,320,000 0.758 
L.P.G. *3.10 3,855,600 1,240,000 0.806 
Refinery Gas Ry 600 


* Figures té eh snes the Oil and Gas Journal for the purposes of this calle elnalats. 


FUEL Rec QUIREMENTS FOR A TOTAL OF F 14, 800 HORSEP POWER Gas T URBINES 


Per Year Per Day Per Hour 

_ (8760 hrs. ) (24 hrs. ) 
Total BTU's required 1.945 x 10 | 5.328x 10 | 2.22x 108 
Barrels of Diesel Oil 328 ,000 898 37.4 
Standard Cubic Feet of Refinery Gas 1. 215: x 10° _5.33 x 108 _139, 000 


With a Thermal Efficiency of 17% each horsepower w ill require 2546 =15,000 BTU's/hour 
<A 
Therefore, Heat Recovered by Steam Generated in Waste Heat Boilers at 46.5% efficiency, 
=.4605 x 15,000 =6,960 BTU’s per hour. 
Note: (a) 1 Horsepower Hour is equivalent to 2,546 BTU’s. 
(b) The Electric Power cost at Sarnia is 0.76 cents per K.W. hour. 





TABLE 9 
DELIVE RY IN MonTHs AS sl st 1956 


Manufacturer 


Gas Turbines operating ee = 
g 


on Diesel oil, or gas In North America In Europe 


\ B Cc D 
1 - 8,500 HP Gas Turbine 12 to 14 . 18 5 to 20 ah itil 12 ey 30° wi 
1 — 6,000 HP Gas Turbine 12 to 14 18 to 20 | 12 36 
1 1,800 HP Gas Turbine | 12to 14 No offering | No offering | No offering 





STEAM TURBINES operating on E F 
600 p.s.i.g. and 125 p.s.i.g. | 


steam 3,200 HP 12 to 14 10 to 11 9 to 10 | 12 to 13 


Although the above table shows little to eae between usted deliveries, the tendency 
would be to place more faith in those quoted for steam than for gas turbines. Centrifugal 
compressor suppliers quote better deliveries than any of the drivers quoted in three cases 
out of four. 


the availability is sure, whereas 
the cost of steam power is apt 
to vary considerably from time 
to time as the cost of refinery 
fuel varies. 


availability are uncertain, their eval- 
uation becomes an art rather than a 
science. 

Accordingly, the design of both 
systems and machines should be flex- 
ible as possible so that varying condi- 
tions and costs may be used to the 
best economic advantage. 


7. Recommendation 

If gas can be made available in the 
quantities required, it is recommended 
that gas turbine drivers be installed. 

This paper shows one of the meth- 
ods used to determine the overall 
economics of a particular problem. 

As in all problems, the solution is 
no better than the factors on which it 
depends. 

When such factors as fuel cost and 
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A Study of Black Liquor Oxidation 


in Towers Packed with 
Asbestos Cement Plates’ 


F. E. MURRAY? 


The process of black-liquor oxidation is treated theoretically and 
an equation is derived for design purposes. Capacity coefficients were 


determined experimentally for two different arrangements of asbestos- 
cement packing sheets. Data include those obtained from a study of a 
full-scale commercial oxidation unit as well as those from an investi- 
gation using a pilot-scale experimental tower with a modified arrange- 
ment of asbestos-cement packing. Results from the pilot-scale study 
indicate that the modified packing yields higher capacity coefficients 
than the packing in the full-scale unit although the mass-transfer co- 
efficients are very similar for both packings using the pilot-scale tower 
under rigorously controlled flow conditions, it was found that the capacity 
coefficients increase with increase in liquor flow rate and are virtually 
unaffected by changes in air-flow rate when this is sufficiently high to 
provide adequate oxygen for the reaction. The significance of the 


experimental results is discussed in detail. 


To oxidation of black cooking 
liquor has been discussed by a 
number of authors (23,4) and is 
now established as a method for odor 
control and improved chemical re- 
covery in kraft pulping operations. 
Although several empirical ap- 
proaches (?, 4) have been used in cal- 
culating the required packing dimen- 
sions for commercial oxidation units, 
no detailed analysis of the combined 
absorption and chemical reaction 
process has been attempted. For the 
evaluation of various packing ma- 
terials and the design of better com- 
mercial towers, it is of considerable 
importance to develop equations for 
the calculation of reliable capacity 
coefficients. 

The initial studies of black-liquor 
oxidation on asbestos-cement plates 
were done on a laboratory scale by 
Wright @), He found that the rate 
of oxidation was very dependent 
upon the pattern on the surfaces of 
the plates and that, at constant air 
flow, over the range of concentra- 
tion studied, the rate was directly 


proportional to the bulk concentra- 
tion of sulphide ion in the black 
liquor. The following mathematical 
treatment yields results in agreement 
with W right’s findings and leads to 
an equation with which capacity co- 
efficients of various tower packings 
may be calculated from experimental 
observations. 

Since oxygen is only very- slightly 
soluble in a hot salt solution, it ap- 
pears reasonable to expect that the 
overall rate of black-liquor oxidation 
will be controlled by the rate of mass 
transfer of oxygen in the liquid 
phase. Accordingly, it is assumed in 
the theoretical treatment that the 
chemical reaction is very rapid com- 
pared to the rates of mass transfer of 
oxygen and sulphide ion to the re- 
action region. 

In the oxidation of black liquor, 
it is highly probable that one chemi- 
cal reaction predominates which may 
be written 


S- + nO2 — Product........ (1) 


Some uncertainty exists regarding 


1Manuscript received December 13, 1957. 
2B.C. Research Council, Vancouver, B.C. 


Contribution from the British Columbia Research Council, 


University of British Columbia, 


Vancouver 8, B.C. This article is based on a paper presented at the 40th Annual Conference 
of The Chemical Institute of Canada, Vancouver, B.C., June 3-5, 1957. 
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the value of n, but for the present 
analysis, it is only necessary that n 
remains constant. According to equa- 
tion (1), the mass of sulphide ion 
oxidized is related to the mass of 

oxygen absorbed by the equation 
Ns = Be 2 ccceadilll 

n 

Using Equation (2) and following 
the methods of Davis and Crandall) 
as outlined by Sherwood and Pig- 
ford (* P- $21), one obtains the equa- 


tion 
De es De HP 
Do nD, . 


ease ipnegls (3) 
H i XL 
ik” Do 





N; = 


An approximation for the term 
D,HP,/nD, can be calculated for 
the absorption of oxygen at 175°F. 
It is reasonable to assume _ that 
D,/D, = 1 and for the oxidation of 
sulphide ion to thiosu!phate, n = 1. 
Based on data available for water (), 
HP, ~ 3 X 10° lb oxygen/Ib water. 
The concentration C ranges from 
about 4 X 10-° to 5 X 10- lb sulphite 
ion/Ib black liquor. For the absorp- 
tion of oxygen in black liquor at 
175°F, the term D,HP,/nD, would 
be somewhat less than for water as 
a result of the lower solubility of 
gases in salt solutions. Consequently, 
this term in Equation (3) is neglig- 
ible for all practical purposes and 
Equation (3) may be written 


H Ss Xi ; 
. 


For fixed gas and liquor flow con- 
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ditions at any given temperature, the 
coefficient of C is substantially con- 
stant and Equation (4) may _ be 
simplified to the form 


Ns = KC. anther eee 


Equation (5) indicates that the 
rate of oxidation of sulphide ion is 
directly proportional to the bulk 
concentration of sulphide ion in the 
solution. This result was demon- 
strated experimentally by Wright ©) 
in his laboratory study. 

In one square foot of packed cross 
section of height dh, there are a 
dh ft? of effective interfacial contact 
area. In this volume of packing, the 
mass of sulphide ion oxidized in unit 
time will be 


N, adh = KC adh........(@) 


The mass of sulphide ion removed 
from the bulk of the solution and 
oxidized in the reaction zone in this 
volume element will be L dC lb/hr. 
Therefore, 


LdC = 


Integration of Equation (7) over 
the full tower height yields i re- 
lationship 


et ee (7) 


Ka = log 


Ka is a capacity coefficient for the 
tower packing expressed in terms of 
sulphide ion oxidized. 

If the oxidation efficiency e is de- 
fined as the percentage of the in- 
coming sulphide ion which is oxid- 
ized in the tower, Equation (8) may 
be written 


Ka = 2303L 100 “ 
_— h . 100 =e * 
Equation (8) or (9) can be used 
to calculate the capacity coefficient 
for any given tower packing when 
h and L are known, and C, and C, 
determined. With Ka knows n, Equa- 
tion (9) can be used to calculate the 
tower dimensions for any desired 
oxidation efficiency. 


Experimental 

Apparatus and Procedures: The 
commercial oxidation unit for which 
data are reported is assembled from 
four sections, each of 4 ft. square 
cross section. The total packed 
height is 34 ft. The packing consists 
of vertical asbestos-cement sheets on 
one inch centres and is arranged to 
give continuous flow surfaces in 
parallel vertical planes from the top 
to the bottom of the tower. A dis- 
tributor at the tower top flows films 
of liquor onto both surfaces of each 
plate. The surfaces are imprinted 
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1 — Pilot-seale 


Figure 
tower packed with modified form of 
asbestos-cement packing. 


experimental 


with a diamond pattern to create 
turbulence in the liquor film as it 
moves downward. In this full-scale 
unit, liquor flow rates were measured 
with a commercial orifice flow meter 
and air flow rates were estimated 
using a velometer. The velometer 
readings had to be averaged because 
of fluctuations in the air flow rate. 
Figure 1 shows the _pilot-scale 
tower which was designed for ex- 
perimental studies on various pack- 
ing types. In this work, the tower 
was packed with corrugated asbestos- 
cement plates which ‘had a screen 
imprint on both major surfaces. The 
packed volume in the tower was ap- 
proximately 1 ft. square in cross- 
sectional area and 7 ft. high. The in- 
dividual sheets were 0.125 in. thick 
and were spaced on 0.66-in. centres. 
The top edge of each sheet was cut 
so that equal amounts of liquid 
would fall on the two sides of the 
sheet when a uniform distribution of 
liquor was directed vertically down- 
ward onto the top of the ‘packing. 
Liquor flow rates were measured 
with a calibrated orifice flowmeter 
in the inlet black-liquor line. Air 
flow rates through the pilot tower 
were regulated by a damper control 
and measured in a 6-in. diameter 
inlet duct with a velometer. A uni- 
form velocity profile in this duct 
made it possible to obtain accurate 
values of G. 

In the study of the full-scale unit, 
normal operating conditions _ pre- 
vailed with liquor flow rates main- 





tained by manual adjustment of 
valves in the inlet feed lines. 


The pilot-scale unit was arranged 
to draw liquor from the inlet line 
of a full-scale oxidation unit and to 
discharge the liquor from its sump 
to the sump of the full-scale towers, 


After changing of either air or 
liquor flow in any experiment, a 
period of 30 minutes was allowed for 
the liquor in the tower sump to reach 
steady-state concentration. Experi- 
ments were randomized ‘) to mini- 
mize the overall effect of slight 
concentration changes which occur 
slowly in the normal operation of 
the mill. 


All chemical analyses were done 
by a colorimetric method “) which 
is the most reliable one available for 
the analysis of sulphide in black 
liquor. 


Results 

Table 1 shows values of the capa- 
city coefficient Ka at various values 
of L in the full-scale oxidation unit. 
The average value of Ka obtained 
from these results is 221 lb/hr ft® 
(unit C). The variance (’) of the ex- 
perimental values is 531 which is 
very close to the known experi- 
mental error variance. There is no 
significant trend in Ka values in 
Table 1 as a function of liquor flow 
rate. 


To obtain reliable information on 
the variation of Ka with liquor and 
air flow rates, a 2 X 3 X 5 factorial 
experiment (8) in duplicate was car- 
ried out using the pilot-scale unit. 
Results of this work are given in 
Table 2. An analysis of variance ‘) 
done on these results shows that 
there is a highly significant increase 
in capacity coefficients with increase 
in liquor flow rate on the corrugated 
packing. The relation between Ka 
and L is given by the equation 
Ka = 39 L°.23, Changing of air flow 
rates over a wide range showed no 
significant effect on Ka values al- 
though at very low air flow rates, 
oxidation rates decrease due to 
oxygen starvation in the tower. This 
behaviour in general has been re- 
ported “), based on studies using a 
small laboratory tower packed with 
raschig rings. 


Discussion 

As was anticipated from a prior 
knowledge of oxygen-absorption in 
general, the rate of oxidation of black 
liquor is controlled by conditions 
in the liquid film and virtually un- 
affected by changes in the gas film. 
This experimental knowledge sup- 
ports the choice of reaction and ab- 
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sorption mechanism in which the 
overall reaction rate is governed by 
the sates of mass transfer of the re- 
actiig components in the liquid 
phase. 

An equation relating K to the 
mass -transfer coefficient k,, for 
oxygen can be readily obtained. 
From Equation (5), N, = KC and 


N r 
from Equation (2), N, = as N, 


can also be expressed in terms of 
k,, by the relationship (° P- *8) 


oo @ Be (Car = Oe ccc ss (10) 


From Equations (10), (2) and (5) 
k, and K are related by the equation 


Since the rate of absorption has been 
shown experimentally to be con- 
trolled by liquid-film conditions, C,; 
can be replaced by HP,. Equation 
(11) then becomes 


Cc 
kL =n (ae) et 


In the oxidation of black liquor, a 
large excess of air is drawn over the 
liquid surface so that P, in the gas 
stream changes very little through- 
out the tower. On the other hand, 
the object of the oxidation is to 
change C from C, to about C,/10 
as the liquor flows from the tower 
top to the sump. Consequently, re- 
ferring to Equation (12), n, H, P, 
and K are essentially constant and 
k,, varies with C. 

As a result of this variation with 
C, k, is not a suitable parameter 
for evaluating packing efficiency in 
black- -liquor oxidation towers. This 
situation is quite different from that 
in processes such as the absorption 
of SO, in an alkaline solution, where 
the concentration of the alkali in the 
solution changes very little whereas 
the partial pressure of SO, in the gas 
stream undergoes a large “change. In 
such a case, the process may be 
treated simply as physical absorption 
without serious error and_ useful 
values of k,, or kg are obtained. 

Despite the unsuitability of k, for 
design purposes in black- “liquor Oxi- 
dation equipment, it is useful to 
calculate approximate average values 
of k,a and H,, for comparative pur- 
poses. From Equation (12), k, = 
10°K (since C/HP, is about 10%). 
From data in Table 2, at L = 9,500, 
Ka = 333 and k,a ~ 333 X 108. 
Based on this value, H,, for the ab- 
sorption of oxygen in black liquor 
on the corrugated asbestos-cement 
packing is about 0.03 ft which, as 
expected, is very much smaller than 


TABLE 1 


CAPACITY COEFFICIENTS AT VARIOUS BLACK LIQUOR FLOW RATES IN A 
FULL-SCALE OXIDATION TOWER 


G = 434 lb./hr. ft.? 


Temp. = 180°F. 


L 2900 | 3312 | 3720 | 4152 | 4968 | 5808 | 6216 | 6624 
Ib./hr. ft.? 
Ka 242 206 226 178 211 233 257 216 
lb. sulphide ion/hr. ft.? 
(unit C) 


TABLE 2 


AVERAGE VALUES OF Ka OBTAINED AT VARIOUS COMBINATIONS OF G AND L. 


Temp. = 175°F. 





G 
Ilb./hr. ft.? 
— 121 322 562 730 845 Averages 
L, 
Ib./hr. ft.? 
honed: 3 i: | S s i 
CT) Se otal 300 293 280 265 273 282 
ee -wigentt aa 4 ae 318 | 318 284 303 
9500... sae Wie aneraa if 348 333 332 315 339 333 
Averages... Ar 321 302 310 299 299 306 
30 
_ 
-_ 
ine 
Ww 
- 4 
20 
10 
4000 6000 8000 


Figure 2—Variation of H.T.U. 
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L (Ib/br.ft? ) 


with liquid flow rates in the oxidation of black 
liquor at 175°F. 
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values of H,, for the purely physical 
absorption of oxygen in water on 
various commercial packings (° P + 200). 


For the process of black-liquor 
oxidation, it is possible to obtain a 
quantity H, equivalent to the height 
of a transfer unit, in terms of mass- 
transfer of sulphide rather than 
oxygen. Values of H, calculated 
from the equation H, = L/Ka are 
shown in Figure 2. It is, however, 
more convenient, at the present time, 
to use Equation (9) for design pur- 
poses. 

A comparison of the thickness of 
the sulphide diffusion film X,, and 
the oxygen liquid diffusion film X;, 
is interesting. Since K = D,/X,, an 
estimate of X,, can be obtained if D, 
is known. D, has a value of about 
7 X 10° ft?/hr at 25°C and can be 
assumed to vary directly as the ab- 
solute temperature and inversely as 
the viscosity of the solution; at 80 ic 
D, ~ 2 X 10+ ft?/hr. Using this 
value of D, and an average value of 
K = 10 (since a on the corrugated 
packing is about 30), X, 0.001 ft. 

Using Equation (12) and the re- 
x and K = Kr Jat 
is found that the ratio X,, is given 
by the equation 


lationships k,, = 


X‘'L/X"_ = = a : (13) 
nC D, 

It was shown previously that HP,/nC 
is approximately equal to 10-*; also, 
D,/D, can be reasonably assumed 
equal to unity. The relationship 
X’,, = 10° X”,, is obtained, so that 
X’,, ~ 10-® ft. or about 10-> in. Since 
X’,, is the distance from the reaction 
zone to the liquid-gas interface, it 
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can be seen that the reaction occurs 
very close to the surface of the 
liquid. 


The convential film-theory was 


used throughout in this work to de- 
scribe the mechanism although pene- 
tration theory ‘'°) may yield a more 
accurate description of the physical 
processes involved. It does appear, in 
any case, that the zone in which the 
chemical reaction occurs is very near 
to the gas-liquid interface. 


Nomenclature 

a = effective surface contact area per 
unit volume, ft.?/ft.? 

C = concentration of sulphide ion in the 


main body of the liquor stream, Ib. 
sulphide ion/Ib. black liquor 

C.i = concentration of oxygen in the 
black liquor at the interface, Ib. 
oxygen/Ib. black liquor 

C; = concentration of sulphide ion in the 
black liquor entering the tower, lb. 
sulphide ion/Ib. black liquor 

Cs = concentration of sulphide ion in the 
black liquor leaving the tower, lb. 
sulphide ion/Ib. black liquor 


D,. = diffusion coefficient for oxygen in 
the solution, consistent units 

D; = diffusion coefficient for sulphide ion 
in the solution, consistent units 

e = percent oxidation efficiency, dim- 
ensionless 

G = air flow rate, lb./hr. ft.? 

h = total height of tower packing, ft. 

H = Henry’s Law constant for solution 


of oxygen in black liquor at the 
prevailing temperature, consistent 
units 

Hy = height of a transfer unit for the 
mass transfer of oxygen in the 
liquid film, ft. 

H, = height of a transfer unit for the 
mass transfer of sulphide ion in 
the liquid film, ft. 


Ka = capacity coefficient for the tower 
packing, lb. sulphide ion/hr. ft.* 
(unit C) 

IK = mass transfer coefficient for sul- 


phide ion in the liquid Ib. sulphide 
ion/hr. ft.? (unit C) 


k, = mass-transfer coefficient for oxy- 
gen in the gas film, Ib./hr. {t.? 
atmosphere 

ki = mass-transfer coefficient for oxy- 
gen in the liquid film, lb. oxygen/ 
hr. ft.2 (unit C) 

L == flow rate of black liquor, Ib./hr. ft.? 

n = number of pounds of oxygen re- 
quired to oxidize one pound of 
sulphide ion, dimensionless 


No = mass of oxygen absorbed in unit 
time over unit interfacial area, 
lb./hr. ft.? 

Ns = mass of sulphide ion transferred to 


the reaction zone in unit time over 
unit interfacial area, lb./hr. ft.? 


P, = partial pressure of oxygen in the 
gas stream, atmospheres 

Xi. = total thickness of liquid diffusion 
films, ft. 

X'L = thickness of the oxygen diffusion 
layer in the liquid, ft. 

X”_ = thickness of the sulphide-ion dif- 


fusion layer in the liquid, ft. 
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Studies on Fermentation Aeration’ 


I. The Oxygen Transfer Coefhcient 


R. E. CARPANE and J. M. ROXBURGH®* 


The transfer coefficients for oxygen from air to water and dilute 
salt solutions and from air to sodium sulphide solutions catalyzed by 
copper ion have been determined in an unbaffled stirred vessel. Up 
to the level of agitation at which the vortex broke the surface there is 
no apparent simple relation between these coefficients. Scale up on the 
basis of sulphite oxidation rates is therefore likely to be successful only 
under conditions of similar geometry and Reynold’ s number. Apparent 
activation energies of 1.2 and 10.6 keal./g.mole were obtained for the 


transfer of 
respectively. 


oxygen 


HE rate at which oxygen is util- 
io by an aerobic fermentation 
is commonly supposed to be limited 
by either or both of the resistances 
to mass transfer encountered, one at 
the gas-liquid interface and one at 
the liquid-cell surface. These two re- 
sistances are in series and the transfer 
rates across them are therefore equal 
under steady state conditions. It is 
probable that each organism has 
characteristic critical oxygen  re- 
quirements and that at utilization 
rates below this critical value growth 
of the organism or product forma- 
tion or both is inhibited. 

The transfer of oxygen from air 
to water or to an aqueous solution 
in which it is sparingly soluble is 
generally accepted to be a “liquid 
film” controlled process“: * 8), The 
transfer coefficient is therefore de- 
pendent on the physical and chemi- 
cal conditions in the liquid phase and 
is not affected by changes in the 
gaseous phase except in so far as they 
affect the solubility of oxygen at 
equilibrium. The rate equation for 
oxygen absorption is then 


to water and 


M/20 sodium sulphite solution 


if the assumption is made that no gas 
phase resistance to mass_ transfer 
exists. The rate of oxygen utilization 
is represented by 


dm 


= K,A,(c)..........-(2) 


indicating that it is dependent, for 
a given organism, on the oxygen 
concentration in the medium, the cell 
concentration (which determines the 
specific interfacial area, A,) and the 
transfer coefficient K,. Both of these 
equations apply only over a limited 
range of dissolved oxygen concentra- 
tion. At higher concentrations cell 
population ‘only becomes limiting, 
characterized by logarithmic — in- 
crease, and at lower concentrations 
growth is inhibited or stopped. 


Bartholomew and co-workers (+ 5) 
and Hixson and Gaden “®) used this 
two stage model in the analysis of 
the effect of aeration on fermenta- 
tions. The methods used resulted in 
values for the product of gas-liquid 
transfer coefficient and_ interfacial 
area (K,A) and the latter workers 
suggested a correlation between this 
product and both power input and 
the superficial air velocity. Quite re- 
cently Johnson and co-workers () 
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Prairie Regional Laboratory, 


have investigated the hydrogenation 
of a-methylstyrene at the surface of 
a_ solid palladium- -alumina catalyst 
using the same two-resistance model. 
They have analyzed their results in 
terms of dimensionless ratios and 
conclude that the transfer coefficient 
at the bubble surface is proportional 
to a power of the gas flow number 
(V.u/o) and a power of a modified 
Reynold’s number (D?Np/u). In ad- 
dition they have concluded that the 
coefficient at the solid surface is de- 
pendent on Reynold’s number alone. 
Their results may be interpreted to 
mean that for a specific fermentation 
in a given fermentor 


KitA « V,9-75 N167,. 0001... (3) 
and 
As tt DR od ceca (4) 


in the terminology of the present 
paper. Equation 3 agrees in general 
with some of the values reported in 
the literature’) and equation 4 
closely approximates the exponent 
of stirrer speed reported by Hixson 
and Baum for the dissolution of 
solids (*) 

Their investigation of the catalytic 
hydrogenation of a-methylstyrene 
shows experimentally that increased 
gas flow rate and increased agitation 
are not equivalent factors in their 
effects on the overall transfer rate. 
Bartholomew and_ co-workers (+ 5) 
made this distinction also in their 
analysis of transfer coefficients in 
streptomy cin and penicillin fermen- 
tations. They suggested that the most 
economic operation is obtained when 
the power input and air flow rate 


are set at minimum values to give 
adequate transfer coefficients at both 
the bubble surface and the cell sur- 
face. 
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Rushton (°, 14) has pointed out the 
significance of the exponent for the 
Rey nold’s number in the analysis of 
gas- liquid contacting. If the expon- 
ent is greater than 0.75, power input 
per unit volume will be decreased 
on scale up, but if the exponent is 
less than 0.75 power imput will be 
increased. For a specific medium in 
a given fermentor Reynold’s number 
is a linear function of impeller speed 
and thus in terms of the notation 
used here 


KitA « N* 


“ 


If the exponent “x” is greater than 
0.75 power imput required per unit 
volume will be decreased on scale 
up. Johnson et al.() report a value 
of 1.67 for this exponent. Cooper, 
Fernstrom and Miller (8), and more 
recently Friedman and Lightfoot (*), 
correlated power imput with aera- 
tion, obtaining exponents of 0.95 and 
unity respectively. 

It is clear that if the two-resistance 
model is accepted, no single factor 
can be considered as an independent 
criterion of the aeration regime in 
a fermentation. For example, in the 
common problem of scale up of 
fermentation, it is not sufficient to 
specify dissolved oxygen concentra- 
tion, c, unless conditions are such 
that the transfer coefficients (K,) 
are comparable, and the cell popula- 
tions (which fix A,) follow the same 
pattern in both cases. 

Direct control of the oxygen up- 
take rate as suggested by Shu either 
by comparison with a “master” fer- 
mentation ‘'S) or by comparison with 
a prev iously determined program (74) 
is one solution to this problem. The 
overall transfer rate is controlled, 
without regard to the individual fac- 
tors, by adjusting the stirring rate. 

Other methods of estimating the 
aeration level in fermentations de- 
pend on measurements of dissolved 
oxygen concentration 4), the respira- 
tion rate of a test fermentation “), or 
the rate of oxidation of a sodium 
sulphite solution catalyzed by a trace 
of cupric ion (> ?* 1), These methods 
rovide an estimate of the product, 
K,A, in either relative or absolute 
terms, and of the effect of operating 
variables on this product. 

Basically, then, the methods of 
measuring aeration for purposes of 
contro] in fermentations fall in three 
categories: 

1) Direct measurement of the rate 
of oxygen utilization 


2) Measurement of the dissolved 
oxygen concentration 


3) Measurement of oxygen trans- 
fer coefficients K,A, (a) directly, 
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from the respiration rate and dis- 
solved oxygen concentration during 
a test fermentation under a variety 
of conditions of operation; (b) in- 
directly, from the relative rates of 
oxidation of sulphite solutions under 
a variety of operating conditions; or 
(c) from a correlation of power 
consumption and _ superficial air 
velocity with K,A, determined as in 
3(a) or (b) above. The first two 
methods make it possible to control 
a fermentation continuously and in 
a fermentor of unknown character- 
istics. Methods in the third category 
depend on a previous knowledge of 
the effect of variables on the char- 
acteristics of the fermentor and do 
not take into account any unusual 
condition encountered during the 
run. On the other hand, the relative 
simplicity of these methods makes 
them attractive if sufficient knowl- 
edge of their vagaries can be ob- 
tained to allow them to be applied 
successfully. Since they are based on 
measurements of K, A it is here that 
information is required. 

The transfer coefficient and inter- 
facial area are affected by stirrer 
speed and air flow rate. Preliminary 
experiments indicated that the com- 
position of the liquid-phase also had 
a marked effect. The present paper 
reports some of the variations in the 
transfer coefficient measured under 
conditions such that the interfacial 
area was known. 


A few values of the transfer co- 
efficient, K,, for the system air-water 
(or oxygen water) ‘have been re- 
ported previously. Davis and Cran- 
dall*) give values of from 0.55 to 
2.4 xX 10-* moles per sq. cm. per 
minute per mole per litre, increasing 
with increasing rate of shear at the 
bubble surface. Schultz and Gaden (?) 
have measured the transfer coeffi- 
cient for the solution of oxygen from 
the air into catalyzed sodium sulphite 
solutions. They obtained values of 
K,, from 9 to 12 X 10-4, but in this 
case lower values were observed as 
the shear rate at the surface in- 
creased. They attributed the unex- 
pected direction of the variation to a 
reduction in the steady state con- 
centration of an unidentified active 
catalyst complex, brought about by 
increased exposure to air. The pro- 
cess was shown to be liquid film 
controlled. 


Experimental 

A thermostatted Pyrex cylinder 
with a flat bladed impeller and no 
baffles was used in all experiments. 
Diameter of the cylinder was 21 cm. 
and of the impeller 17.8 cm. Four 
to six litres of solution were used in 


each run. This is roughly the equi- 
valent of the — used by 
Schultz and Gaden (), Sulphite solu- 
tions were in all cases 0.05 M Na,SO, 
catalyzed with CuSO, (0.001 M). 
The pH was adjusted to 8.0 before 
oxidation. Samples were added to an 
excess of iodine and the excess back- 
titrated with standard sodium thio- 
sulphate (N/10) to a starch end 
point. At stirrer speeds high enough 
to cause substantial depression of the 
surface of the liquid at the center, 
an area correction was made by as- 
suming the surface to be conical. 
The correction was essentially neg- 
ligible under all but the most extreme 
conditions. 


Dissolved oxygen was determined 
by the W inkler method recom- 
mended in “Standard Methods for 
the Examination of Water and Sew- 
age” (17), Saturation values were de- 
termined in all cases and compared 
with values given in this reference 
and in the literature. 


A number of samples were taken 
during each run and _ statistical 
methods used, where applicable, 
obtain the best fitting line from 
which K,, was then calculated. In 
calculating K,, for sulphite solutions 
the assumption was made that the 
concentration of dissolved oxygen 
was zero at all times. Equation 1 then 
reduces to 


K dm/dt (5) 

a a finite ae eee ) 
A.C 

The saturation value for oxygen in 

the sulphite solution was assumed to 

be the same as for water at the same 

temperature and pressure. 

The transfer coefficient for water 
and salt solutions was calculated 
from the integral form of the rate 
equation 

1 Ge = Gis 


Ki = “a + aeO) 
- A (te = t) A co = Cee . 


either graphically or numerically. 


Results and Discussion 

The transfer coefficients for 
oxygen to water, 0.5 M_ sodium 
chloride and 0.05 M sodium sulphite 
are shown in Figure 1 as a function 
of stirrer speed. The coefficients for 
water and for sodium chloride solu- 
tion are “normal” in that they in- 
crease more or less regularly. The 
exponent at low r.p.m. for water is 
about 0.55 and for the sodium chlor- 
ide solution 0.75. These increase 
sharply at just below 100 r.p.m. to 
about 2.6. It is interesting to note 
that O’Connor and Dobbins ') using 
an oscillating stirrer found that the 
transfer coefficient varied as_ the 
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squere root of the frequency of 
osciliation both for water and sul- 
hite solutions. This correlation 
agree d with their theoretical analysis 
of the rate of surface replacement 
which was predicted to be propor- 
tional to the square root of the mean 
velocity fluctuations. The inflection 
point in Figure 1 corresponds rough- 
ly to the point at which the vortex 
becomes appreciable. Rushton (") has 
suggested that the transfer coefficient 
is a function of Reynold’s number, 
D*Np 


’ 


in baffled tanks or at low 
shear rate in unbaffled ones, but that 
. DN? 
the Froude number, N , must also 
g 


be considered under conditions lead- 
ing to vortex formation. 


The significant feature here is not 
the exact correlation with stirrer 
speed or the dimensionless numbers 
describing the correlation, since the 
reaction vessel bears little resem- 
blance to a fermentor. Rather the 
important feature is the very rapid 
rise in the value of the transfer co- 
efficient as agitation is increased. It 
can be inferred from these values for 
K,, that a coefficient of at least 10-° 
moles/cm.2. min. (mole/litre) and 
perhaps much higher is to be ex- 
pected in a fermentor. The correla- 
tion for K,;, proposed by Johnson 
and co-workers (), with the expon- 
ent of Reynold’s number 1.67 for the 
hydrogen — a-methylstyrene system 
may be applicable in a fully baffled 
tank, but further experimental work 
is required to confirm this inference. 


The effect of stirrer speed on the 
transfer coefficient for oxygen to the 
sulphite solutions shown in Figure 1 
is quite different. Qualitative agree- 
ment with the work of Schultz and 
Gaden (2) is indicated but their values 
for K,, are nearly twice as large. The 
decrease in value from 7.6 X 10-* to 
36 X I0* as agitation is increased 
can only be caused by a decrease in 
the rate of the controlling step in 
the complex oxidation reaction for 
sulphite. This may very well be a 
reduced concentration of a catalyst 
intermediate as suggested by Schultz 
and Gaden. Polarographic analyses 
of 0.001 M copper sulphate solutions 
with and without sulphite present 
show respectively two electron and 
one electron oxidations with a slight 
shift in the half wave potential. This 
suggests that the copper ion is pre- 
sent as a complex with sulphite as 
might be expected. The species re- 
sponsible for this polarographic wave 
is not the active catalyst complex 
since the height of the wave increases 
in proportion to the concentration 


Slope 0.55—=> 
1.0 


0.5 





0 30 40 





Water —>/ 
a 


if 


an Na Cl 


£ OY Slope 0.75 








60 80 100 120 160 200 


RPM 


Figure 1—The effect of stirrer speed in an unbaffled vessel on the transfer 
coefficient for oxygen from air to water and to sodium chloride and sodium 
sulphite solutions. 


of copper ion added but catalytic 
activity does not. Since the value of 
the transfer coefficient at the highest 
stirrer speed drops below the value 
that would be expected for diffusion 
into a salt solution without reaction, 
the assumption of zero oxygen con- 
centration in the bulk of the solution 
may not be valid. If this is so, then 
all of the values of K,, for sulphite 
solutions are low. 

In Figure 2 the effect of additional 
solute on the transfer coefficient for 


0.05 M_ sodium sulphite solutions 
without agitation is shown. The 
small disturbances from convection 


currents and sampling do not affect 
the results, and in fact the stirrer can 
be operated up to about 25 r.p.m. 
with only slight variation in the 
measured transfer coefficients. Sod- 
ium acetate and glucose depress the 
transfer rate considerably, sodium 
chloride has about the same effect 
as in water alone and no depression 
is observed with ethanol or addi- 
tional sodium sulphite. Since these 
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values are not corrected for the 
solubility of oxygen in the solution 
the transfer coefficients may actually 
be higher than shown. Apparently 
the assumption that the solubility of 
oxygen is the same in 0.05 M sodium 
sulphite as in water is satisfactory 
since increasing the concentration to 
0.20 M (5 gm./100 ml.) did not alter 
the transfer coefficient calculated 
using this assumption. Similar values 
were obtained with agitation but at 
the higher solute concentrations 
foaming caused erratic results. The 
effect of additives on transfer of 
gases to water has been investigated 
before “° 2°) mainly with reference 
to surface active agents, with rather 
variable results. Cullen and David- 
son (1°) indicate that pure sodium 
dodecyl sulphate has no effect on the 
rate of transfer of carbon dioxide to 
water but the addition of a little 
dodecyl alcohol to the dodecyl sul- 
phate depresses the rate markedly. 

The effect of temperature on the 
transfer coefficients is shown in 
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Figure 2—The effect of additional solutes on the transfer coefficient for oxygen 
from air to a M/20 sodium sulphate solution. No agitation. 
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Figure 3—Arhenius plot of the effect of temperature on the transfer coefficient 
for oxygen from air to water and to a M/20 sodium sulphate solution. Stirrer 
speed 30 r.p.m. 
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Figure 3. The apparent activation 
energies calculated from the slopes 
are 10.3 kcal./g. mole for 0.05 M sul- 
phite and 1.2 kcal./g. mole for water, 
In the latter case the temperature co- 
efficient should be closely related to 
changes in the diffusion coefficient 
with temperature. There is both 
theoretical (+, 21,22) and experimen- 
tal (#8, 24) evidence that the transfer 
coefficient varies as the square root 
of the diffusion coefficient. The 
effect of temperature on this dif- 
fusion coefficient is apparently not 
too well established, but four values 
reported (25, 26, 27,28) predict activa- 
tion energies for the transfer coeffi- 
cient of from 1.2 to 3.3 kcal./g.moles. 
The Stokes-Einstein equation pre- 
dicts a square root dependence of the 
transfer coefficient on viscosity and 
inversely on the square root of ab- 
solute temperature. Other empirical 
relations for the transfer coefficient 
are similar and all predict activation 
energies somewhat higher than that 
reported here, but of the same order 
of magnitude. 


Conclusions 


The transfer coefficients calculated 
from sulphite oxidation rates bear no 
simple relation to the coefficients for 
the transfer of oxygen from air to 
water under the same conditions. 
From the limited range of variables 
covered in this investigation the fol- 
lowing general conclusions can be 
drawn: 

1) Transfer rates to water and 
non-reactive salt solutions increase 
rapidly with increased stirring rates. 
In the unbaffled vessel used the co- 
efficient increases as the 2.6 power 
of stirring rate above the point at 
which a vortex begins to form. 

2) Transfer rates to a catalyzed 
sulphite solution decrease with in- 
creased stirring rate. For a still solu- 
tion the coefticient is about fifteen 
times the value measured for water 
decreasing to one-half the value for 
water at the highest stirring rates 
tested. 

3) The effect of temperature on 
the transfer coefficient is much great- 
er for sulphite solutions than for 
water. Apparent activation energies 
are 10.6 and 1.2 kcal./g.mole respec- 
tively. 

If sulphite oxidation rates are to 
be used successfully in the control 
and scale up of aeration in fermenta- 
tions it is quite evident that com- 
parison must be made at the same 
Reynold’s number. Only then is it 
likely that the relation between 
transfer coefficients for the test sul- 
phite solution and for the fermenta- 
tion medium ‘will be the same. UI- 
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timately, this requires similar geom- 
etry and fully baffled conditions. 
Further, there is no reason to sup- 
pose that a linear relation exists 
between sulphite oxidation rates and 
the effective aeration level for a 
fermentation run under the same 
conditions. Sulphite oxidation rates 
reflect changes in aeration in only 
a qualitative way unless the nature 
of the relation between them can be 
established. 


Nomenclature 


dm 
—- = rate of oxygen transfer moles/ 
dt min. 


A = interfacial area, air-liquid, cm.?/ 
cm.$ 

As = interfacial area, solid-liquid,: 
cm.?/cm.3 

x = a constant 

c = concentration of oxygen in the 
bulk of the liquid, moles/litre 

Cs = saturation concentration of oxy- 


gen, moles/litre 
Ct1, Cte = concentration of oxygen at times 


ti, te 
Ki = oxygen transfer coefficient, gas- 
ar F : mole 
liquid, moles/cm.2 min. {| —— 
litre 
Ks = oxygen transfer coefficient, solid- 
is j : mole 
liquid, moles/cm.? min. {| — 
litre 





g = gravitational constant, cm./min.? 
D = impeller diameter, cm. 
p = density, gm./cm.? 
be = viscosity, gm./cm./min. 
o = surface tension, gm./min.? 
N = agitator speed, r.p.m. 
Vs = superficial air velocity, cm./min. 
D?Np ; ‘ 
= Reynold’s number, dimension- 
b less 
DN? ’ ? 
= Froude number, dimensionless 
g 
Vee : ; 
—- = gas flow number, dimensionless 
o 
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The Preparation of a Crystalline 
High Explosive of Controlled Particle Size 
by Precipitation with Water 


from Acetone Solution’ 


A. M. PENNIE? 


The preparation of fine RDX of controlled particle size has been 
studied and a satisfactory continuous process has been developed. One 
of the main features of the process is the use of an impinging jet mixing 
head which obviates many of the variables encountered in previous 


methods of preparation. 


RYSTALLINE high explosives in a 

finely divided form are impor- 
tant ingredients in a variety of pro- 
pellant and detonator compositions. 
The size limits of the crystals required 
are critical and it is essential for satis- 
factory performance of the finished 
compositions that there should be 
good batch to batch reproducibi. ity of 
the crystalline ingredients. Wet mil- 
ling of the crystals in a disk mill fails to 
give the desired size and precipitation 
from a solvent offers the most attrac- 
tive alternative. 

The particular high explosive con- 
sidered for the above applications was 
RDX*. Previously the method of 
preparation of very fine RDX to meet 
propellant and detonator requirements 
was never too satisfactory. Generally, 
fine RDX was prepared by the r rapid 
drowning of hot RDX-acetone solu- 
tion in excess water with stirring. 
When manufactured in this manner 
the product always varied over a wide 
range—2 to 100 uw. If very fine material 
was desired, say, 2-8 u, then the entire 
product had to be laboriously sieved 
and the over-size rejected. Sieving of 
such fine material is far from satis- 
factory and if carried out by hand 
shows. great variation from operator 
to operator. During this rapid drown- 
ing method of precipitation conditions 
are far from constant and it is not 
surprising to find the product ranging 
over wide limits. The concentration 
of acetone increases as the drowning 
progresses and the rate of precipita- 


*Cyclotrimethylenetrinitrainine. 


tion of RDX by water decreases as 
the water to acetone ratio decreases. 

The type of stirring employed has 
always been arbitrary and depended a 
great deal on local conditions and on 
the plant available. It is always diffi- 
cult to reproduce stirring character- 
istics and this contributed greatly to 
the differences in procedure and 
results associated with fine RDX pro- 
duction. The term “good agitation” 
is rather a nebulous one which has 
been used freely to cover this rather 
important factor. 

A review of the general method of 
preparation of fine RDX by rapid 
drowning of a hot charge of acetone- 
RDX into cold water indicates that it 
is influenced by several factors: 


(a) The rate of addition of RDX 
solution varies as the level of 
the liquid in the header tank 
falls. 

(b) The concentration of the 
water-acetone mixture in the 
precipitator changes  contin- 
uously during the agitation. 

(c) The temperature in the pre- 
cipitator increases as the RDX 
solution is added. 

(d) Stirring characteristics change 
as the volume of the mixture in 
the precipitator increases. 

All these factors can be brought 

under control if a system of contin- 
uous precipitation is adopted. 


Continuous precipitation 
A continuous pot precipitator was 
assembled and_ several exploratory 


1Manuscript received February 14, 1958. 


2Chief Superintendent, Suffield Experimental Station, Ralston, Alta. 
Note: This article is based on work carried out at the Canadian Armament Research and De velop- 
ment Establishment, Valcartier, Que., and is covered by Canadian Patent 533487 and British Patent 


741756. 
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runs were carried out to obtain infor- 
mation as to how the spread of crystal 
size was affected by the temperature 
of feed, the temperature of mixing, 
speed of stirring and method of addi- 
tion of the feeds. A 10% solution of 
RDX-acetone and water was metered 
into the vessel which was equipped 

with cooling coils and twin propeller 
stirrers and charged initially with the 
equilibrium mixture of aqueous ace- 
tone. Samples of the product were 
withdrawn and examined under a 
biological microscope, and the crystal 
size determined by measurement with 
an eyepiece micrometer at a magnifi- 
cation of 430. Results from the many 
runs performed indicated that a high 
degree of reproducibility from run to 
run was possible. Although the range 
of crystal size was controllable, the 
spread within the r range was unaccept- 
able. 

Three methods of addition of the 
RDX-acetone and the diluent water 
were employed: 

(a) Water and RDX-acetone dia- 
metrically opposite on the sur- 
face. 

(b) Water at the bottom of the 
pot and RDX-acetone on the 
surface. 


(c) Water and RDX-acetone at 
varying depths below the sur- 
face (as in Figure 1). 


Results 


Tables 1 and 2 show that with a 
fixed feed rate, increasing the stirrer 
speed results in an increase in fines 
and raising the temperature of the 
precipitating water and the reactor 
increases the crystal size. 

Table 3 indicates that increasing the 
stirrer speed has practically no effect 
on the crystal size since the majority 
of mixing and precipitation has taken 
place by virtue of the two streams 
merging at the point of addition. 
However, increasing the temperature 
of the reactors again results in an 
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incrcase in crystal size but the depth 
of te addition point below the sur- 
face has little or no effect on the 
partic cle size. At 30°C this method of 
mixing always gives a product rang- 
ing from 20-80u irrespective of the 
stirrer speed. This is indicative that 
dilution and precipitation are more 
rapid with this method of addition. 

Ex perience with the continuous 
precipitator indicates that the spread 
of crystal size arises from the poor 
mixing of the streams of the two 
components. With this type of mix- 
ing, before equilibrium is reached, 
acetone rich and water rich layers or 
pockets are formed in the body of 
the liquid with intermediate concen- 
trations between the two extremes. 
These layers may only exist for a 
very short time but the time is suffi- 
cient for precipitation and crystal 
growth to take place. The acetone 
rich layers will naturally give crystals 
of a larger size than the water rich 
layers and the intermediate layers will 
produce crystals of sizes varying 
between the two extremes. 

This problem of layers or pockets 
of varying concentrations can be 
eliminated ‘considerably by having the 
mixing take place in a large number 
of places at the same time. Equi- 
librium concentration is thus more 
rapidly reached. The layers, under 
such conditions, will then be smaller, 
the time of precipitation will be 
shorter, and hence the time of crystal 
growth will be decreased. All this 
should tend to give small crystals of 
approximately the same size. 


Precipitation by spray mixing 
Apparatus 

As a result of the observations in 
the continuous pot precipitator it was 
decided to emp'oy a system of 
multiple mixing points. Accordingly, 
a spray head consisting of a series of 
impinging jets was constructed. The 
apparatus is illustrated in Figures 2 
and 3. The body was made of alum- 
inum and the jets of brass threaded 
into the body. The two sets of jets— 
one for the RDX-acetone solution 
and the other for the diluent water, 
are mounted at 90° to each other. The 
outer duct, carrying the RDX-ace- 
tone solution, is surrounded by a hot 
water jacket. All pipes carrying the 
RDX-acetone solution were hot water 
jacketed. 

It was hoped that this system of 
impinging jets would make it possible 
to control the crystal size by adjus- 
ting the rates of flow of the i imping- 
ing streams. If the resultant streams 
are made acetone rich the crystals 
should be large and if the resultant 


WATER 


RDX -ACETONE 


<i 





Figure 1—Method of addition of RDX-Acetone and water at varying depths 
below the surface of the reactor. 


TABLE 1 
EFFECT OF METHOD OF ADDITION ON PARTICLE SIZE 
Addition of both components (diametrically opposite) on the surface of the liquid 


RDX-Acetone Temp. of Pptn. water Temp. of Temp. of RPM = Crystal 
feed RDX-Acetone | gm./min. pptn./water pptor. _ stirrer Size 

gm./min. "€ re a bb 
1140 : 54 1140 9 21 400 60-100 
1140 54 1140 9 21 600 20-100 
1140 54 1140 9 21 800 5-100 

| | 
1140 54 | 1140 30 30 400 | 80-180 
1140 54 1140 30 30 600 40-140 
1140 54 1140 30 30 800 | 20-120 
TABLE 2 


EFFECT OF METHOD OF ADDITION OF PARTICLE SIZE 
Addition of water at the bottom of the pot and RDX-Acetone on the top surface 


RDX-Acetone| Temp. of Pptn. water Temp. of | Temp. of RPM | Crystal 


feed RDX-Acetone | gm./min. pptn./water  pptor. stirrer size 
gm. ‘min. ss we "e be 
1140 | 54 1140 9 21 400 | 30-100 
1140 54 1140 9 21 | 600 20-100 
1140 54 1140 9 21 | 800 5-100 
| | 

1140 51 | 1140 = | 27 | 31 400 | 40-160 
1140 51 1140 31 31 600 | 20-140 
1140 | 52 | 140 | 31 31 | 800 | 20-100 
1140 | 54 1140 9 21 600 20-100 
1140 53 1140 32 30 600 10-120 
1140 53 1140 44 40 600 10-180 
10-200 


1140 53 1140 55 53 600 
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Figure 2—Section through the mixing head showing the solution and water 
reservoirs, the heating jacket and the spray nozzles. 
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JACKET WATER OUTLET 


DETACHABLE NOZZLES 


Figure 3—Plan view of the spray mixing head showing the detachable nozzles. 


RDX-ACETONE WATER 


e Si fle 


CONTROL ne CONTROL VALVE 


FLOW METER FLOW METER 
COLLECTING FUNNEL 
SPRAY MIXING HEAD 

Figure 4—Flow sheet illustrating the process. 


—-WATER 


RDX ~ ACETONE 





/2” DIA 


Figure 5—Apparatus employed when using jets of 1% in. diameter and greater. 


streams are water rich the crystals 
should be smaller. 


Operation of the process 


The flow sheet (Figure 4) indicates 
the main features of the system. Hot 
filtered ten percent RDX-acetone was 
fed through a calibrated flowmeter 
via the jacketed pipe lines to the 
mixer head. Precipitating water at the 
required temperature was likewise fed 
to the mixer through a calibrated 
flowmeter. Ten percent RDX-acetone 
solution was chosen as a convenient 
concentration to handle since it 
would not cause any serious crystal- 
lization troubles in the lines should 
the temperature fall slightly below 
the operating limit. In an experimental 
unit such as this, stoppages for inspec- 
tion, measurements and adjustments 
were frequent, and if a more concen- 
trated solution had been used the risk 
of “icing-up” would have increased. 

The spray was allowed to run for 
several minutes so that equilibrium 
could be achieved before a sample 
was taken. Generally speaking, the 
runs for sampling purposes were 
short, five to ten minutes, but in some 
cases runs of forty to fifty minutes 
duration were performed and _ the 
results agreed with those of shorter 
duration. By varying the flow of both 
the feeds effective control over a wide 
range of the mass ratio of the RDX- 
acetone solution to water was ob- 
tained. Jets of different sizes were 
employed and the effect of varying 
the nozzle diameter was examined. 
The stream from the mixer head was 
directed into a funnel and run to a 
centrifuge for filtration. Earlier in 
this investigation a vacuum bag filter 
was used, but with such a fine pro- 
duct filtration was very slow. As a 
result of this holdup on the filter 
there was a tendency for the crystals 
to grow in the supernatant liquor. 
Use of the centrifuge speeded up 
filtration and obviated this short- 
coming. 


Results 


Table 4 indicates the range of 
optimum conditions for the produc- 
tion of very fine uniform crystals. 
This lies between mass ratios RDX- 
acetone/water of 0.1:1-—30.3:1. It will 
be noted that the temperature of the 
water used as precipitant is not 
critical, except at temperatures some- 
what above normal room temperature. 
Under the latter conditions, the rate 
of precipitation is slower and there is 
an opportunity for crystals to con- 
tinue growing, thus giving a wider 
range of particle size. In addition it 
can be seen that for mass ratios of 
RDX-acetone to water between 0.4:1, 
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and 3:1 a product of larger crystals 
over wider ranges of particle size is 
obtained under similar operating con- 
ditions as the mass ratio is increased. 

From Table 5 it is seen that for low 
mass ratio of RDX-acetone to water 
the effect of increasing the nozzle 
size is not particularly significant. 
Very fine crystals (1-2 ) are formed 
using nozzles with holes from 1/16” 
diameter to 5/32” diameter under 
similar conditions of mass ratio and 
temperature. However, at mass ratios 
greater than 1:1, larger crystals and 
wider ranges of ‘particle size are ob- 
tained as the size of the nozzle hole 
is increased. 

To investigate the effect of nozzles 
larger than 5/32” diameter an ar- 
rangement as shown in Figure 5 was 
adopted. 

Experiments with this equipment 
using nozzles of 4” diameter indicated 
that although the finest crystals were 
not obtained it was still possible to 
prepare fine (2-10 mw) crystals with 
regularity and reproducibility. As in 
the case of the multiple-spray head 
work the crystal size could be con- 
trolled by adjusting the mass ratio of 
the feeds and the temperature of the 
precipitating water. 


Conclusions 

The system of impinging jets can 
be used to prepare very fine crystals 
of closely control’ed particle size. 
The method employed is independent 
of the many variables previously 
encountered in other mixing processes 
and in addition to being continuous it 
eliminates the indefinable quantity of 
“stirring”. 

No sieving is required and the 
process can be operated with the 
minimum of labour. The resultant 
mother liquor for recovery is reason- 
ably strong and concentration of ace- 
tone for reuse is simple and satisfac- 
tory. 

With the smaller units, the mixer 
head is capable of producing about 10 
pounds per hour of the fine (1-2 w) 
RDX; with a pair of the 4” nozzles 
an output of 240 pounds per hour is 
possible. Clearly a battery of mixer 
heads could be operated efficiently to 
obtain any desired production capa- 
city, and the output from both these 
units could be increased by using a 
14% solution of RDX in acetone. 
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TABLE 3 
EFFECT OF METHOD OF 


ADDITION ON PARTICLE SIZE 


Addition of both components below surface 


Temp. of Pptn. water Temp. of 


RDX-Acetone ; gm./min. pptn./ water 
oe °C 
24% below surface 
52 1140 31 
52 1140 30 
52 1140 30 
52 1140 40 
314” below surface 
52 1140 30 
52 1140 31 
52 1140 30 
52 1140 41 


TABLE 4 


Nozzle Temp. of | Total feed 
diameter Water gm./min. 
inches "] 
1/16 13 3100 
1/16 15 2300 
1/16 4 2700 
1/16 10 2700 
1/16 15 2700 
1/16 32 3240 
1/16 50 2700 
1/16 15 2700 
1/16 4 1820 
1/16 13 3380 
1/16 15 2825 
1/16 10 1200 
1/16 14 2850 
1/16 13 3200 
1/16 14 3400 
1/16 15 3000 
1/16 15 4000 
| 1/16 15 4000 
1/16 13 3600 
1/16 15 3450 


percent solution. 


TABLE 5 


Temp. of 


pptor. 


32 
30 
30 
41 


RPM 


stirrer 


400 
600 
800 
600 


400 
600 
800 
600 


WATER 0.1:1 — 3:1 





Crystal 
size 
bu 
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20-80 
40-100 


20-90 


20-80 
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40-100 


Crystal Size 
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COMPARISON OF Mass RATIO AND NOZZLE SIZE 





Nozzle Temp. of | Total feed 
diameter Water gms./min. 
inches “Cc 

1/16 10 2700 
3/32 9 2700 
| 5/32 10 | 2700 
1/16 13 | 3200 
3/32 | 9 2400 
5/32 10 2400 
1/16 is | 4000 
3/32 | 9 4000 
5/32 | 10 | 4000 
1/16 | 13 | 3600 
3/32 | 9 | 3600 
5/32 10 3600 
1/16 15 3450 
3/32 9 3240 
5/32 10 3240 
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The Phlegmatization 


A. M. PENNIE? and T. 


of Fine RDX’ 


A successful method for phlegmatization of fine RDX has been 
developed. Small quantities for the phlegmatizing wax are dissolved 
with RDX in acetone and the wax is co-precipitated as a very thin film 


on individual RDX erystals. 


This method is incorporated into the technique for precipitation of 


fine RDX previously reported. 


HE method of preparing RDX* 
‘bon rstals described in a previous 
paper ‘!) produces material of fine 
particle size w hich on drying, consti- 
tutes a possibie hazard due to inherent 
sensitiveness and a tendency to dust- 
ing particularly when it has to be 
handled in the dry state. Accordingly, 
consideration was given to dev eloping 
a method of phlegmatizing or desensi- 
tizing this fine material. 

A general method by which crys- 
talline high explosives are phlegma- 
tized with wax is to slurry the 
explosive and wax with water at a 
temperature above the melting point 
of wax. With violent agitation a 
single dispersed phase of explosive 
crystals enrobed in globules of 
molten wax is formed. On cooling 
the slurry, the globules retain their 
identity as pellets of solid wax, each 
enclosing one or more crystals of 
explosive, and can be readily sepa- 
rated from the water dispersant. The 
pellets may vary in size from several 
tens of microns to about one inch in 
diameter depending on such factors 
as the size of the explosive crystals, 
the relative concentration of explosive 
to water in the slurry, the type of 
vessel used and the type of agitation 
employed. 


®Cyclotrimethylenetrinitramine. 


Similar methods may be used to 
coat an explosive with viscous densi- 
tizing oils, 

The conventional method described 
was not considered suitable for use 
with fine RDX since it was desired to 
achieve the phlegmatization with as 
little phlegmatizer as possible in order 
to retain the maximum explosive 
potential and to ensure that the ph'eg- 
matized material was free flowing and 
discrete—not agglomerated with oil 
or wax. 

Co-precipitation of the wax and 
RDX by water from an acetone solu- 
tion was considered as a_ possible 
method. It was thought that the wax, 
in. precipitating, might form on the 
surface of each RDX crystal rather 
than as an inclusion within the crystal 
or as discrete wax particles. 


Experimental 

Two waxes were investigated — 
Beeswax and “Acrawax”* both of 
which are soluble in hot acetone to 
the extent of one percent or more. In 
all cases the RDX employed was 
RDX(B)** of World War II manu- 


facture. In this investigation the 


*Glyco Products. 

*°RDX(B), manufactured by the Bachmann 
Process, contains 6-12 percent of HMX 
(cyclotetramethylenetetranitramine ). 


1Manuscript received February 14, 1958. 


2Chief Superintendent, Suffield Experimental Station, Ralston, Alta. 

3Canadian Armament Research and Development Establishment, Valcartier, P.Q. 
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- STERLING* 


weight of wax used varied from 5 to 
10% of the dry weight of RDX; ice. 
from 0.5 to 1% of the weight of 
RDX-acetone solution when a 10% 
solution of RDX-acetone was em- 
ployed. 

The acetone was heated to about 
50°C in a dissolving kettle and the 
required amount of wax added in a 
finely divided state in order to facili- 
tate solution. When all the wax was 
dissolved, the dry RDX was added 
and the temperature maintained at 
50°-54°C to complete solution. The 
technique of co- precipitation by 
spraying with impinging jets was as 
described in a previous paper () and 
the size of the nozzles was varied 
from 1/16 to 4 in. diameter. The 
resultant slurry was run directly to 
a centrifuge, washed and subsequently 
dried in a drying stove. A sample run 
is recorded in Table 1. 


Results 
Ratio-solution/precipitant 
The results shown in Table 2 indi- 
cate that for a given solution/precipi- 
tant ratio and a fixed nozzle size, the 
product was similar to that obtained 


from runs in which no wax was 
added. 


Impact sensitivity 

One of the important tests applied 
to explosives is the determination of 
their sensitivity to impact —i.e. the 
distance in centimcters through which 
a falling weight must drop on a 
sample of the material to cause it to 
explode. For the fa-ling weight or 
impact test a 2 kiogram weight is 
generally employed, and the result, 
the impact sensitivity, is represented 
in kilogram centimeters. When com- 
parisons against a known standard are 
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Acetone 98% 


TABLE 1 


OPERATING CONDITIONS AND Data OBTAINED FROM A TyPICAL RUN 
UsING » RanOwAR AS PHLEGMATIZING AGENT 


100 parts 


RDX 11 parts 


Beeswax 

Solution Spray Jets 
Water Spray Jets 
Solution Feed Rate 


Water Feed Rate 


0.7 parts (6.4% of weight of RDX) 
6 jets each of 3/32 in. diameter 

6 jets each of 3/32 in. diameter 
2.65 lbs./min. 


5.30 Ibs./min. 


Temperature of RDX-Acetone Solution 52°C 


Temperature of Water 10°C 


Resultant Particle Size of RDX: 
Specific Surface 


Impact Sensitivity 


2— Su 
23,000 cm.?/cm.* 


300% of F. of I. 
— Su. 


71 kg.cm. for 2kg.wt. 
for uncoated fine RDX 2 





“% wax on crystals So 


TABLE 2 
EFFECT OF VARYING THE RATIO-SOLUTION/ PRECIPITANT 


Range of Crystals 








made the term Figure of Insensitivity 
is employed. 





Trial Ratio ee ee San SE ks ae 
No. Solution/Precipitant | 
_| Waxed Unwaxed 
mn Me 
1 0.5 :1 2-5 2-5 
2 0.57:1 2-5 2-7 
3 0.67:1 2-5 2-10 
4 0.8 :1 5-10 5-12 
5 a; oo 5-15 5-15 
6 2 10-40 10-50 
TABLE 3 
IMPACT SENSITIVITY RESULTS 
Crystal \ Impact % of F. of I. 
Samples | Size Wax | Wax Used Sensitivity | rel. to standard 
Mv | kg. cm. 
Standard 2-8 0 | none 24 | 100 
Coarse 5 
Unwaxed >150 0 none 17 70 
Waxed >150 10 Beeswax 51 | 212 
Fine No. 1 2-5 5 Beeswax 71 296 
2 2-5 8.5 Beeswax 51 212 
$ 2-5 8.5 Beeswax 48 200 
4 2-10 5 Beeswax | 58 242 
5 2-10 8 Acrawax 48 200 
6 2-10 5 Acrawax | 45.5 189 
7 15-40 0 none | 14 61 
8 15-40 5 Beeswax | 48 200 








shown in Table 3 are calculated from 
50% detonations with a 2 kg weight. 


Impact sensitivity tests were applied The RDX used as a standard was 


from the same lot employed in the 


to the product after drying. These 
tests were carried out on the Im- 
proved Bruceton Machine. The values 


wax precipitations. 
Two coarse samples 


(>150u), 
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waxed and unwaxed, are included in 
the table for purposes of comparison. 

Reference to Samples 7 and 8 shows 
that while RDX/wax of 15-40u has a 
Figure of Insensitivity (F. of I.) of 
200% of the 2-8u standard, its F. of I. 
relative to unwaxed RDX of 15-40u 
size is over 300%. 

It can be seen that coating fine 
RDX with 5% of wax results in a 
marked decrease in sensitivity, and 
that the use of more than 5% of wax 
is of little advantage. The results may 
be summarized by stating that for 
material 2-—510u, sensitivity is de- 
creased by a factor of 2. For material 
15-+40u, the sensitivity is decreased 
by a factor of 3. 


Effect of various waxes 

Two waxes were employed—Bees- 
wax and Acrawax. The results re- 
corded in Table 3 indicate that there 
is no appreciable difference in the 
phlegmatizing properties of the two 
and it is thought that any compatible 
wax soluble in acetone and capable of 
wetting a clean RDX surface, could 
be employed successfully. 


Static charge on fine RDX 

Fine RDX which has not been 
ph'egmatized by wax does not readily 
pick up a static charge. It was ob- 
served, however, that fine RDX 
which had been phlegmatized with 
beeswax as previously outlined tended 
to become statically charged when 
handled after drying. Fine RDX 
phlegmatized with Acrawax did not 
become statically charged. 


Filtering problem with fine RDX 

Unphlegmatized fine RDX is read- 
ily separated from the acetone-water 
solution after precipitation by means 
of a perforate basket centrifuge. 
When RDX of the same crystal size is 
ph'egmatized by co- precipitation with 
wax, filtration becomes quite a 
problem. Not all the wax is precipi- 
tated onto the RDX, and the small 
amount left after precipitation is 
suspended as a very fine froth and 
tends to clog the filter medium of the 
centrifuge. Use of an imperforate 
basket centrifuge greatly improved 
the separation. 


Microscopic examination of 
phlegmatized RDX 

(a) Under White Light 

There is no apparent difference 
between phlegmatized and unphleg- 
matized RDX of the same crystal 
size when examined under white light. 
Reflection from the faces of coated 
and uncoated crystals appeared sim- 
ilar. No discrete pellets of wax could 
be observed. 
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(b) Under Ultra Violet Light 

It was impossible to observe any 
difference in individual crystals, 
coated or uncoated, but when a small 
aggregate of fine phlegmatized mater- 
ial was examined and compared with 
a similar quantity of unwaxed fine 
RDX of the same crystal size, it was 
possible to detect a slight fluorescence 
from the phlegmatized material. 


Flotation test for free wax 

Although no individual wax par- 
ticles could be seen on microscopic 
examination, it was thought that some 
might be present in the product. Both 
waxes have a density less than that of 
water, while RDX has a much higher 
value. To test for the presence of wax 
particles, samples of fine phlegmatized 
RDX were wetted with water con- 
taining a small amount of a surface 
active agent and then stirred into a 
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large volume of water. In no case did 
wax particles appear to separate from 
the RDX. 


Analysis for non occluded wax 

Free partictes of wax and surface 
coated wax were determined by 
means of extraction with carbon 
tetrachloride in which the solubility 
of RDX is negligible but the wax is 
readily soluble. In all cases it was 
found that the wax value obtained 
was about one percent lower than the 
amount of wax used in the spraying 
solution, calculated on the basis of 
dry RDX. 


Conclusions 

A method for producing fine RDX 
phlegmatized with wax has been de- 
veloped. This method is simple in 
operation and readily adaptable to 
use on a plant scale. From the impact 
sensitivities, flotation tests and wax 


analyses, it is evident that the wax, 
which is co-precipitated with the 
RDX, is deposited on the surface of 
the crystal. It has not been possible to 
determine whether all surfaces of a 
crystal receive a coating of wax, al- 
though from the nature of precipita- 
tion, this would be a_ reasonable 
assumption. It should be noted, how- 
ever, that for fine RDX* desensitized 
with 5% of its weight of wax, to coat 
all surfaces would result in an average 
thickness of wax film of about 4 x 
10-® cm. 


°2-—>5UM size has a specific surface of 25,000 
cm2/cm3, 
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